THE ASTROPHYSICAL JOURNAL 


AN INTERNATIONAL REVIEW OF SPECTROSCOPY AND 
ASTRONOMICAL PHYSICS 


VOLUME 86 OCTOBER 1937 NUMBER 3 


ON THE MASSES OF NEBULAE AND OF 
CLUSTERS OF NEBULAE 


F. ZWICKY 


ABSTRACT 


Present estimates of the masses of nebulae are based on observations of the /umi- 
nosities and internal rotations of nebulae. It is shown that both these methods are 
unreliable; that from the observed luminosities of extragalactic systems only lower 
limits for the values of their masses can be obtained (sec. i), and that from internal 
rotations alone no determination of the masses of nebulae is possible (sec. ii). The 
observed internal motions of nebulae can be understood on the basis of a simple me- 
chanical model, some properties of which are discussed. The essential feature is a central 
core whose internal viscosity due to the gravitational interactions of its component 
masses is so high as to cause it to rotate like a solid body. 

In sections iii, iv, and v three new methods for the determination of nebular masses 
are discussed, each of which makes use of a different fundamental principle of physics. 

Method iii is based on the virial theorem of classical mechanics. The application of 


this theorem to the Coma cluster leads to a minimum value M =4.5X 10'°M@ for the 
average mass of its member nebulae. 

Method iv calls for the observation among nebulae of certain gravitational lens 
effects. 

Section v gives a generalization of the principles of ordinary statistical mechanics to 
the whole system of nebulae, which suggests a new and powerful method which ulti- 
mately should enable us to determine the masses of all types of nebulae. This method 
is very flexible and is capable of many modes of application. It is proposed, in par- 
ticular, to investigate the distribution of nebulae in individual great clusters. 

As a first step toward the realization of the proposed program, the Coma cluster of 
nebulae was photographed with the new 18-inch Schmidt telescope on Mount Palomar. 
Counts of nebulae brighter than about m= 16.7 given in section vi lead to the gratifying 
result that the distribution of nebulae in the Coma cluster is very similar to the dis- 
tribution of luminosity in globular nebulae, which, according to Hubble’s investiga- 
tions, coincides closely with the theoretically determined distribution of matter in 
isothermal gravitational gas spheres. The high central condensation of the Coma 
cluster, the very gradual decrease of the number of nebulae per unit volume at great 
distances from its center, and the hitherto unexpected enormous extension of this 
cluster become here apparent for the first time. These results also suggest that the 
current classification of nebulae into relatively few cluster nebulae and a majority of 
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field nebulae may be fundamentally inadequate. From the preliminary counts reported 
here it would rather follow that practically a// nebulae must be thought of as being 
grouped in clusters—a result which is in accord with the theoretical considerations of 
section v. 


In conclusion, a comparison of the relative merits of the three new methods for the 
determination of nebular masses is made. It is also pointed out that an extensive 
investigation of great clusters of nebulae will furnish us with decisive information re- 
garding the question whether physical conditions in the known parts of the universe 
are merely fluctuating around a stationary state or whether they are continually and 
systematically changing. 

The determination of the masses of extragalactic nebulae con- 
stitutes at present one of the major problems in astrophysics. Masses 
of nebulae until recently were estimated either from the luminosities 
of nebulae or from their internal rotations. In this paper it will be 
shown that both these methods of determining nebular masses are 
unreliable. In addition, three new possible methods will be out- 
lined. 


I. MASSES FROM LUMINOSITIES OF NEBULAE 


The observed absolute luminosity of any stellar system is an in- 
dication of the approximate amount of luminous matter in such a 
system. In order to derive trustworthy values of the masses of 
nebulae from their absolute luminosities, however, detailed informa- 
tion on the following three points is necessary. 

1. According to the mass-luminosity relation, the conversion fac- 
tor from absolute luminosity to mass is different for different types 
of stars. The same holds true for any kind of luminous matter. In 
order to determine the conversion factor for a nebula as a whole, 
we must know, therefore, in what proportions all the possible lumi- 
nous components are represented in this nebula. 

2. We must know how much dark matter is incorporated in nebu- 
lae in the form of cool and cold stars, macroscopic and microscopic 
solid bodies, and gases. 

3. Finally, we must know to what extent the apparent luminosity 
of a given nebula is diminished by the internal absorption of radia- 
tion because of the presence of dark matter. 

Data are meager’ on point 1. Accurate information on points 2 


«It should, however, be mentioned that certain spiral nebulae seem to be stellar 
systems similar in composition to the local Kapteyn system of our galaxy. For such 
systems the conversion factors may with some confidence be set equal to the con- 
version factor of the Kapteyn system. See also E. Hubble, Ap. J., 64, 148, 1920. 
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and 3 is almost entirely lacking. Estimates of the masses of nebulae 
from their observed luminosities are therefore incomplete and can 
at best furnish only the lowest limits for the values of these masses. 


II. MASSES FROM INTERNAL ROTATIONS OF NEBULAE 


It has apparently been taken for granted by some astronomers? 
that from observations on the internal rotations good values for the 
masses My of nebulae could be derived. Values of the order of 
My = 10? Mo up to My = 4 X 10” Mo were obtained in this 
way,’ where Mo = 2 X 10% gr is the mass of the sun. A closer 
scrutiny of the behavior of suitably chosen mechanical models of 
stellar systems, unfortunately, soon reveals the fact that the masses 
of such systems, for a given distribution of average angular velocities 
throughout the system, are highly indeterminate, and vice versa. 
This conclusion may, for instance, be derived from the consideration 
of two limiting models of a nebula as a mechanical system. 


A. MODEL OF A NEBULA WHOSE “INTERNAL VISCOSITY”’ IS NEGLIGIBLE 


This model consists of a heavy and small nucleus of mass M, 
around which a given number, 7, of stars of average mass M, K 
M,/n describe planetary orbits. The mutual gravitational interac- 
tions between these outlying stars are negligible, and the system 
may therefore be said to have an internal “‘viscosity”’ equal to zero. 
It is obvious that under these circumstances we may build up models 
that satisfy almost any specifications in regard to total mass, total 
luminosity, and internal distribution of luminosity as well as dis- 
tribution of the average angular velocities. We may, for instance, 
distribute our stars over the six-dimensional manifold of all pos- 
sible planetary orbits (including the epochs or phases) in such fashion 
that the average angular velocity of the resulting system S, is zero 
in every point. Since all these orbits are essentially non-interacting, 
we may reverse the sense of rotation (direction of stellar motion) 
in an arbitrary number of these orbits. In this way a system S of 
specified distribution of average angular velocities may be construct- 
ed whose remaining characteristics, such as the mass, the luminosity, 


2 E. Hubble, The Realm of Nebulae (New Haven: Yale University Press, 1936), p. 
179; also Ap. J., 69, 150, 1929. 


220 F. ZWICKY 


and the external shape, are identical with those of S,. Thus, the 
observed angular velocities in themselves give no clue regarding the 
mass of the system. 


B. MODEL OF A NEBULA WHOSE INTERNAL VISCOSITY IS VERY GREAT 
This model is built up of stars, dust, and gases in such fashion 
that the gravitational interactions between the various components, 
as well as direct impacts, influence the path of every component 
mass in a radical way. Many changes in energy and momentum of 
every component mass will take place during time intervals that are 
short compared with the time an unperturbed mass of the same 
initial velocity would consume to traverse the system. 

Conditions of motion in this model are analogous to the condi- 
tions of motion of elementary particles in a star. This model of a 
nebula, therefore, will rotate like a solid body, regardless of what its 
total mass and the distribution of mass over different regions of the 
system may be. The conclusion which has sometimes been put for- 
ward,” that constant angular velocity necessarily implies uniform 
distribution of mass, is obviously erroneous. Furthermore, it is again 
seen that the rate of rotation of a stellar system has no very direct 
bearing on its total mass. 

C. ACTUAL NEBULAE 

Good mechanical models of actual nebulae may presumably be 
constructed by combining the distinctive features of the two limit- 
ing cases described in the preceding sections. Such a combined mod- 
el will possess a central, highly viscous core whose relative dimen- 
sions are not negligible but are comparable with the extension of the 
whole system. If the outlying, and among themselves little inter- 
acting, components of the nebula had no connection with the central 
core, we might, at a given instant, observe average angular velocities 
Q which, as a function of the distance r from the center of rotation, 
would be given by 


Q(r) = Q, = const. forr <1, (1) 
where 7, is the radius of the core. For r > r,, the angular velocity 


Q(r) would be essentially arbitrary. In reality, however, the vis- 
cosity will not drop abruptly to zero at r = r,. From an inspection 
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of the distribution of the outlying masses in many nebulae it would 
seem that these masses at some previous time must have formed part 
of the central core. They may have been ejected from this core 
because they acquired high kinetic energy through many close en- 
counters, or they may be the result of a partial disruption of 
the core by tidal actions caused by encounters with other nebulae 
(Jeans). At the moment these outlying masses have passed the 
boundary of the core, their average tangential velocities must have 
been of the order r,Q,. Assuming that these masses in regions r > 7, 
are essentially subject only to the gravitational attraction of a cen- 
tral spherical core and that the interactions among themselves may 
be neglected (the internal viscosity of the outlying system is equal 
to zero), the average angular velocity Q(r) outside the core will be 
approximately given as 


Q(r) = forr > fe. (2) 


This relation simply expresses the fact that a mass m which, on be- 
ing ejected from the core with a tangential velocity », relative to this 
core, describes an orbit whose angular momentum, 


mrw = m[rovt + (3) 


isa constant. If the average 7; for many particles leaving the core is 
zero, the average angular velocity @ = & of all the particles in a 
given point is obtained from 


rw = rQ(r) = (4) 


which is the same as (2). Unfortunately, relation (2) is superficially 
very similar to the relation obeyed by the angular velocities w, of a 
system of circular planetary orbits around a heavy central mass M. 
For such an orbit we have 


mre? = (5) 
or 
1/2 
= = (6) 


| 


292 F. ZWICKY 


where I is the universal gravitational constant. The similarity of 
the dependence on r in (2) and (6) will in reality become still greater, 
since the rate of decrease of Q(r) with increasing values of r will in 
actual nebulae be more gradual than that given by (2), owing to the 
fact that the internal viscosity will not vanish abruptly at r = 7, 
but will disappear gradually with increasing r. The observed angular 
velocities in the outlying regions of nebulae actually show a de- 
pendence on r which resembles the relation (6). This relation was, 
therefore, sometimes erroneously used for the determination of the 
mass M. The preceding discussion, however, indicates that the ob- 
served angular velocities may be adequately accounted for on the 
basis of the considerations resulting in relation (2) rather than in 
relation (6). This again shows clearly that it is not possible to de- 
rive the masses of nebulae from observed rotations without the use of 
additional information, since the relation (2) does not contain the 
mass M at all. 
D. FURTHER DISCUSSION OF THE ROTATION OF NEBULAE 


From the analysis of an appropriate mechanical model of actual 
nebulae we have derived in the preceding section the following ap- 
proximate dependence on ¢ of the average tangential velocity V;: 


= forr <1 (7) 
forr > 1%. (8) 


This dependence is pictured in Figure 1 . The observed curvesV ,(r) 
are similar to the schematic curve of Figure 1. From this similarity 
we infer that the cores of nebulae possess considerable internal vis- 
cosity. The problem of deducing theoretical values for this viscosity 
therefore arises. It will be interesting to see whether the gravita- 
tional interactions in dense stellar systems are sufficient to account 
for the fact that such systems rotate like solid bodies, or if it is neces- 
sary to introduce matter in the form of dust particles and gases. 
The point of view adopted in the preceding discussion automati- 
cally eliminates the discrepancy which was thought to exist, with 
respect to the distribution of mass in nebulae as derived from data on 
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internal rotations, on the one side, and from the shape and the dis- 
tribution of the isophotal contours, on the other. The tremendous 
increase of surface brightness from the edge, r = ro, of the core of 
nebulae to their center, 7=0, indicates a correspondingly large in- 
crease of mass density. The erroneous idea? that the constancy of the 
angular velocity throughout the core necessitates the assumption of 
a constant mass density therefore created an apparently insoluble 
paradox. This paradox, however, disappears as soon as we introduce 
the idea of an internal gravitational viscosity of stellar systems, 
which equalizes the angular velocity throughout such systems re- 
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Fic. 1.—Velocity of rotation in nebulae 


gardless of the distribution of mass. How this viscosity may be ex- 
pressed in terms of the gravitational interactions of stars will be dis- 
cussed in another place. 

One further interesting problem presents itself. The distribution 
of matter in a stellar system which has an internal viscosity as high 
as we have assumed it to be in the cores of nebulae should rapidly 
converge toward stationary conditions. The determination of the 
density distribution in such a system should be analogous to the 
determination of the density distribution in gravitating gas spheres. 
R. Emden’s analysis? of such spheres may, therefore, prove useful 
in the study of globular and elliptical nebulae as well as in the study 
of cores of spiral nebulae. 

In concluding this section the question may be raised whether 


3 Gaskugeln (Leipzig, 1907). 
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data on the internal rotation of nebulae, supplemented by certain 
additional information, make possible the determination of nebular 
masses. This question, in principle, must be answered in the affirma- 
tive. For instance, data on internal velocities combined with the 
virial theorem discussed in the next section may ultimately furnish 
good values for the masses of globular and elliptical nebulae as well 
as the cores of certain spirals. In the case of open spirals, an investi- 
gation of the geometrical structure of the spiral arms, combined with 


a 


Fic. 2.—Rotation of nebulae and shape of spiral arms 


velocity data, promises to be helpful. In this connection a few re- 
marks may be in order regarding the possibility of the formation of 
spiral arms in our simple mechanical model of nebulae. 

Consider the two alternatives pictured in Figure 2: (a) a non- 
rotating core, that is, Q, =o, and (6) a rotating core, 2, o. Sup- 
pose that some disturbance, such as tidal action, causes the ejection 
from the core of anumber of stars at P. The resulting orbits of three 
stars, 1, 2, and 3, are schematically indicated in Figure 2. For Q, #0 
distinct spiral arms may be formed, provided a sufficient number of 
stars are ejected whose radial and tangential velocity components 2, 
and », relative to the core are in certain favorable relations to the 
peripheral velocity r,Q, of the core. 

It should be remarked that the orbits 1, 2, 3, etc., in Figure 2b 
are ellipses only if outside the core the corresponding masses are no 
longer acted upon by the forces which originally caused their ejec- 
tion from the core. If, for example, the ejection is brought about by 


| | 
i] 
\ i 
/ ‘\ 
b 
) 
| 
} 


THE MASSES OF NEBULAE 225 


the close encounter of the nebula with another nebula, tidal actions 
will continue until the second nebula has moved far away and our 
simple model must be correspondingly modified. It may suffice here 
to point out that in addition to the creation of angular momentum in 
a nebula, a close encounter may also set up radial pulsation in this 
nebula. This pulsation will not, in general, be spherically symmetri- 
ca]. Asa first approximation we can schematically characterize it by 
what might be called the “tidal ellipsoid” or “ellipsoid of pulsation.” 
The magnitudes and signs of the radial velocities which are induced 
along the principal axes of this ellipsoid by the passing second nebula 
approximately determine the character of the resulting pulsation in 
the nebula under consideration. 

The probable existence of pulsation in many nebulae further com- 
plicates the interpretation of differences in the observed velocities 
V, along the line of sight of various parts of a nebula. If we have 
pulsation, V, differs from the value of V, owing to rotation alone, by 
amounts which depend on the position of the tidal ellipsoid relative 
to the two vectors which define the line of sight and the axis of rota- 
tion. In two obvious, simple cases the functions V,(7) may graphi- 
cally assume the forms represented by the two dotted curves (V, and 
in Figure tr. 

The possible existence of pulsation also suggests that we cannot, 
without additional knowledge, interpret ellipsoidal shapes of nebulae 
as being due to rotation alone. Values of the masses of nebulae 
which are derived on the assumption that elliptical nebulae are 
stellar systems in rotational equilibrium must therefore be viewed 
with suspicion. 

Before embarking on any detailed analysis of spiral structures it 
will be advisable to obtain first qualitative information concerning 
the shape of the spiral arms in relation to the sense of rotation of the 
core. Promising results in this direction have already been secured 
by Dr. W. Baade of the Mount Wilson Observatory, who has found 
that the tips of the spiral arms of the Andromeda nebula are curved 
in the direction determined by the sense of rotation of the core, as 
pictured in Figure 2b. I am indebted to Dr. Baade for a private com- 
munication of this unpublished result. 

Summing up, we may say that present data on internal rotations 
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furnish, at best, minimum values M,,;, for the masses of the nebulae. 
Such minimum values are obtained if we assume that nebulae are 
stable systems whose components have velocities v inferior to the 
velocity of escape v, from the system. Therefore 


Ve > Umax 2 (9) 
But 
1/2 
(10) 


where M,, 7., and {, are the mass, the radius, and the angular veloc- 
ity of the core. Consequently, the mass M of the nebula is 


M > M, > (11) 


For example, in the case of NGC 4594 we have ro > 4.3 X 1077 cm 
and 7,.Q, = 4 X 107 cm/sec, which gives 


M > 5 X 10% gr = 2.5 X 10°Mo. (12) 


The data used are taken from the paper on NGC 4594 by F. G. 
Pease,‘ who, in 1916, measured the spectroscopic rotation of this 
nebula. Pease also was the first to point out that the central parts 
of nebulae rotate like solid bodies. 

Since the determination of nebular masses is of considerable im- 
portance in modern astrophysics, three new methods for the solution 
of this problem are briefly outlined in the following pages. Two of 
these I have proposed already in previous communications.*° The 
third, in a restricted form, has been applied thus far only to clusters 
of stars.’ Its applicability to clusters of nebulae remains to be in- 
vestigated. The purpose of this paper is to discuss the basic princi- 
ples on which the new methods rest. 


4 Proc. Nat. Acad., 2, 517, 1916. 

5 F. Zwicky, Helv. physica acta, 6, 110, 1933. 

°F. Zwicky, Phys. Rev., 51, 290, 1937, and §1, 679, 1937. 
7H. von Zeipel, Jubilaeumsnummer d. A. N., p. 33, 1921. 
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III. THE VIRIAL THEOREM APPLIED TO CLUSTERS OF NEBULAE 


If the total masses of clusters of nebulae were known, the aver- 
age masses of cluster nebulae could immediately be determined from 
counts of nebulae in these clusters, provided internebular material 
is of the same density inside and outside of clusters. 
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Fic. 3.—The Coma cluster of nebulae 


As a first approximation, it is probably legitimate to assume that 
clusters of nebulae such as the Coma cluster (see Fig. 3) are mechan- 
ically stationary systems. With this assumption, the virial theorem 
of classical mechanics gives the total mass of a cluster in terms of the 
average square of the velocities of the individual nebulae which con- 
stitute this clusters But even if we drop the assumption that clus- 
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ters represent stationary configurations, the virial theorem, in con- 
junction with certain additional data, allows us to draw important 
conclusions concerning the masses of nebulae, as will now be shown. 
Suppose that the radius vector from a fixed point in the cluster to 
the nebula (c) of mass M, is 7,. For the fixed point we conveniently 
chose the center of mass of the whole cluster. The fundamental law 
of motion of the nebula (¢) is 
M, =F,, (13) 
where F, is the total force acting on M,. Scalar multiplication of this 
equation with r, gives 


d? > dr.\? 
1 zd = . - 
2 dt? (M or2) 'o F, + u.(@*) . (14) 


Summation over all the nebulae of the cluster leads to 


Vir + 2Kr, (15) 


where 9 = >-M,r?2 is the polar moment of inertia of the cluster, 


. > . . . 
Vir = dor, Fe is the virial of the cluster, and Kz is the sum of the 


kinetic energies of translation of the individual nebulae. If the clus- 
ter under consideration is stationary, its polar moment of inertia 0 
fluctuates around a constant value 9,, such that the time average 
of its derivatives with respect to time is zero. Denoting time aver- 
ages by a bar, we have in this case 


Vir = — 2Kr. (16) 


On the assumption that Newton’s inverse square law accurately de- 
scribes the gravitational interactions among nebulae, it follows that 


Vir = E,, (17) 


E,=—- <y (18) 


where 


o,v lov 


| 
| 
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is the total potential energy of the cluster due to the gravitational 
interactions of its member nebulae. Equation (16) thus takes on the 
well-known form 


E, = 2Kr= Mai = > Mai, (19) 
where v, is the velocity of the mass M,. In order to arrive at a quan- 
titative estimate of the total mass . / of a globular cluster of nebulae, 
we assume as a first approximation that these nebulae are, on the 
average, uniformly distributed inside a sphere of radius R. In this 
case 

E 3r. 
(20) 


We may also write 


Mv = Me, (21) 


where the double bar indicates a double average taken over time and 
over mass. Therefore, from (19), (20), and (21), 


_ 
M = (22) 
This relation (22) can also be derived if we take the time average 
of equation (14) which holds for an individual nebula. The time 
average of the left side of (14) disappears if the mass M, is a member 
of a stationary system. Thus we have 


= — 2ker = — (23) 


where vir, is the virial of the nebula (¢) and &,7 is its kinetic energy. 
Now, in a sphere of mass . # and radius R, the density p, if uniform, 


is equal to p = 3.M@/4nrR3. The force F,(r) which acts on M, is 


therefore 


F,(r) R3 (24) 
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(26) 


Since it has been assumed that all the nebulae combined produce a 
uniform distribution of matter throughout the sphere, the average 
nebula spends equal times in equal volumes, and we have 


X = —, 2 


where the double bar again designates a double average with respect 
to time and mass. Consequently, 


M Br’ (28) 


as before. 

From the distribution of the brighter nebulae in the Coma cluster 
pictured in Figure 3 it is apparent that the assumption of uniform 
distribution is not fulfilled. But it is also evident that the actual 
potential energy £, will have a value which, in order of magnitude, 
is correctly given by (20). Even if we crowded all the cluster nebulae 
into a sphere of radius R/2, the value of E, would only be doubled. 
Also, we can increase E, only very little if we assume that the M,’s 
run through a wide range of values. For instance, if we assumed that 
practically the whole mass . / of the cluster were concentrated in 
only two or three nebulae of mass .#/2 or .@ /3, respectively, and 
that these masses had mutual distances as small as R/10, we should 
arrive at values for the potential energy: 


— — 307.02 


—~ 


which, combined with (23), results in | 
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These values are of the same order of magnitude as in (20). The fol- 
lowing inequalities must therefore be considered as conservative 
estimates for the possible maximum values of the average kinetic 
energy and the minimum values of the total mass . /: 


— 
2aKr = —E,< - (30) 
R 
and 
Re 


We apply this relation to the Coma cluster of nebulae whose radius 
is of the order of 2 X 10° light-years. From the observational data 
we do not know directly the velocities v of the individual nebulae 
relative to the center of mass of the cluster. Only the velocity com- 
ponents 2, along the line of sight from the observer are known from 
the observed spectra of cluster nebulae. For a velocity distribution 
of spherical symmetry, however, we have 


v= 302. (32) 
Therefore 
(33) 


From the observations of the Coma cluster so far available we have, 
approximately,’ 


v2 = 5 X 10'cm? sec™. (34) 


This average has been calculated as an average of the velocity 
squares alone without assigning to them any mass weights, as actual- 
ly should be done according to (21). It seems, however, as Sinclair 
Smith*® has shown for the Virgo cluster, that the velocity dispersion 
for bright nebulae is about the same as that for faint nebulae. As- 
suming this to be true also for the Coma cluster, it follows that the 


J., 83, 499, 1936. 
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mass-weighted means of v? and the straight means are essentially the 
same. Furthermore, in calculating (34) we have used velocities 
which belong to the bright nebulae, since only these have been meas- 
ured. If brightness can be taken as a qualitative indication of mass, 
the error in substituting (34) for (21) cannot be great. We must, 
nevertheless, remember that, strictly speaking, the determination of 
_M dy the virial theorem is subject to the difficulty of calculating 
through the application of an averaging process which involves the 
as yet unknown masses. The mass . % as obtained from the virial 
theorem, can therefore be regarded as correct only in order of mag- 


nitude. 
Combining (33) and (34), we find 


M> 9 X (35) 


The Coma cluster contains about one thousand nebulae. The aver- 
age mass of one of these nebulae is therefore 


M > 9 X 10% gr = 4.5 X 10° Mo. (36) 


Inasmuch as we have introduced at every step of our argument in- 
equalities which tend to depress the final value of the mass . % the 
foregoing value (36) should be considered as the lowest estimate for 
the average mass of nebulae in the Coma cluster. This result is 
somewhat unexpected, in view of the fact that the luminosity of an 
average nebula is equal to that of about 8.5 X 107 suns. According 
to (36), the conversion factor y from luminosity to mass for nebulae 
in the Coma cluster would be of the order 


Y = 500, (37) 


as compared with about y’ = 3 for the local Kapteyn stellar system. 
This discrepancy is so great that a further analysis of the problem 
is in order. Parts of the following discussion were published several 
years ago, when the conclusion expressed in (36) was reached for the 


first time.’ 
We inquire first what happens if the cluster considered is not st-a 
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tionary, in which case the virial theorem (16) must be replaced by 
one of the two inequalities 


2Kr + Vir Z0, (38) 


where the bars denote averages taken over time intervals which are 
comparable with the time it takes one nebula to traverse the whole 
system. The smaller sign needs no further consideration, since in- 
stead of resulting in equation (19) it leads to the inequality 


— E, > 2Kr. (39) 


This, in turn, means that the inequality (33) is further enhanced, 
and we arrive at a lowest value for the mass .# greater even than 


(36). 
The other alternative is 


— E, < 2Kr. (40) 


We may combine this inequality with the principle of conservation 
of total energy E of the whole cluster 


Kr + E, =E, (41) 
which, added to (40), leads to 
E>-—Kr. (42) 


This means that the cluster is expanding. In particular, if E = owe 
have 


E,= — Kr (43) 


instead of E, = —2Kyz, for a stationary system. Equation (43) 
means that the cluster will ultimately just fly apart. In this case the 
mass . # of the cluster still has a value of half that arrived at in (36), 
and the discrepancy between the conversion factors y and 7’ remains 
materially the same in order of magnitude. 

If we wish to reduce the mass . #/still further so as to make ap- 
proximately y = 7’, we must put E >oand Ky > — E,. By as- 
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suming this, however, we run into two serious difficulties. In the 
first place, it is difficult to understand why under these circumstances 
there are any great clusters of nebulae remaining in existence at all, 
since the formation of great clusters by purely geometrical chance is 
vanishingly small, as will be shown in another paper. In the second 
place, the cluster nebulae, after complete dispersion of a cluster 
would still possess velocities practically identical with their original 
velocities. The field nebulae in general, which under the assumed 
circumstances can hardly have a velocity distribution different from 
that of former cluster nebulae, should therefore have a dispersion in 
peculiar velocities comparable to that of cluster nebulae. The aver- 
age range of peculiar velocities among field nebulae, however, seems 
to be of the order of 150 km/sec only. This observation, if correct, 
excludes values of the total energy E which are sufficiently greater 
than zero to reduce the discrepancy between y and y’ to a satisfac- 
tory degree.’ It will, nevertheless, be advisable to obtain more data 
on the velocities of both cluster nebulae and field nebulae in order 
to arrive at accurate values of the dispersion which characterizes the 
respective velocity distribution functions. 

In addition it will be necessary to develop methods which allow 
us to determine the relative amounts of internebular material in 
clusters as well as in the general field. 

It should also be noticed that the virial theorem as applied to clus- 
ters of nebulae provides for a test of the validity of the inverse square 
law of gravitational forces. This is of fundamental interest because of 
the enormous distances which separate the gravitating bodies whose 
motions are investigated. Since clusters of nebulae are the largest 
known aggregations of matter, the study of their mechanical be- 
havior forms the last stepping-stone before we approach the inves- 
tigation of the universe as a whole. 

The result (36) taken at face value of course does not mean that 
the average masses of field nebulae must be as great as those of clus- 
ter nebulae. From the general principles discussed in a following 
section one would rather expect the heaviest nebulae to be favored 
in the process of clustering. 

The distribution of nebulae in the Coma cluster, illustrated in 
Figure 3, rather suggests that stationary conditions prevail in this 
cluster. It is proposed, therefore, to study the Coma cluster in more 
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detail. On the other hand, the virial theorem can hardly be used with 
much confidence in cases such as the Virgo cluster and the Pisces 
cluster.? These clusters are much more open and asymmetrical than 
the Coma cluster and their boundaries are thus far ill defined. Ac- 
curate values of the gravitational potentials in these clusters are diffi- 
cult to determine. 

In passing it should be noted that the mechanical conditions in 
clusters of nebulae are in some important respects different from the 
conditions in clusters of stars. During close encounters of stars only a 
minute part of their translational energy is transformed into internal 
energy of these stars, if the extremely rare cases of actual impacts 
are disregarded. Nebulae act differently. In the first place, close en- 
counters and actual impacts in a cluster of nebulae must occur dur- 
ing time intervals which are not very long compared with the time 
of passage of one nebula through the entire system. Therefore, a 
considerable tendency exists toward equipartition of rotational and 
internal energy of nebulae with their translational energy. Star clus- 
ters in some ways are analogous to gas spheres built up of monatomic 
gases, whereas clusters of nebulae may be likened to gas spheres 
built up of polyatomic gases. This difference in internal character- 
istics may lead ultimately to serious consequences, as will be seen 
from the following line of reasoning. 

We again start from the virial theorem (19), in the form in which 
it may be written for stationary gravitational systems, namely, 


E,+ 2Kr=o. (44) 


Admitting that on close encounters of nebulae kinetic energy of 
translation Ky; may be transformed into energy of rotation Kr and 
internal energy E, of nebulae, we must replace the restricted form 
(41) of the energy principle by the more general equation 


Kr+Kr+Er+£,=E. (45) 
Combining (44) and (45), we have 


= —Kr+Kr+Er. (46) 


9F. Zwicky, Proc. Nat. Acad., May, 1937. 
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Encounters among nebulae tend to establish equipartition among 
translational, rotational, and internal energies of nebulae, analogous 
to the equipartition of energy among different degrees of freedom 
which is so well known in ordinary statistical mechanics. The me- 
chanical conditions in a cluster of nebulae should converge toward 
a state for which 


where A > o. Forinstance, \ = 1 if we include in £; only the kinetic 
energy of the first three fundamental modes of pulsation of the whole 
system. The older the system becomes, the more of its many de- 
grees of freedom may be expected to share in the equipartition of 
energy. Other happenings excluded, the value of \ would, therefore, 
increase almost indefinitely as the average value E; of E; is deter- 
mined for time intervals of increasing length. In any case, admit- 


ting (47), 
E=)Kr>o. (48) 


Weare thus led to the uncomfortable conclusion that the total energy 
of astationary cluster of nebulae should be positive. This is a con- 
tradiction in itself, since it means that on the assumption of sta- 
tionary conditions in a cluster we have proved that the cluster really 
cannot be stationary but must ultimately fly apart. It would seem, 
therefore, that clusters of nebulae analogous to gravitational gas 
spheres which are built up of polyatomic gases could never repre- 
sent stationary configurations. In reality, however, condition (47) 
cannot be reached. In contradistinction to the rotational energy of 
polyatomic molecules at low temperatures, the rotational energy kp 
of a nebula cannot become equal to its observed translational energy 
kr. Long before the equipartition (47) among the various types of 
energy can be established by close encounters among nebulae, these 
nebulae will have been partially or completely disrupted. We are 
here confronted with processes which are analogous to the dissocia- 
tion of polyatomic molecules when their average kinetic energy of 
translation—that is, the temperature of the gas—becomes too high. 

The preceding considerations open up interesting new vistas on 
the change in time of nebular types in clusters where close encoun- 


Kr=Kp. Er= Xz, (47) | 
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ters are much more frequent than among field nebulae. In the first 
place, the most compact and most massive nebulae are presumably 
the least vulnerable to disruption. These nebulae, therefore, are 
destined to survive the longest. This effect of selective elimination of 
nebular types may, in part, be responsible for the difference in the 
representation of types among cluster nebulae as compared with 
field nebulae. In the second place, we should expect a considerable 
number of stars, as well as matter in dispersed form from disrupted 
nebulae, to be scattered through the internebular spaces within clus- 
ters. Sufficiently large amounts of internebular matter in clusters 
might seriously change our estimate (36) of the average value of 
nebular masses as derived from the preceding application of the virial 
theorem to clusters of nebulae. It is therefore the intention to under- 
take a series of observations which may throw some light on the 
problem of the density of internebular matter in clusters, as com- 
pared with the density of matter in the general field. Until such ob- 
servations have been made it will be well to keep in mind that, al- 
though the determination of average nebular masses from the virial 
theorem may be viewed with considerable confidence, this method is 
not entirely free from objections which have not yet been satisfac- 
torily dealt with. 

In principle the virial theorem may also be applied to describe the 
mechanical conditions in an individual nebula. Actually a direct ap- 
plication is difficult, since it is not possible to measure separately, 
as in the case of a cluster of nebulae, the velocities of the individual 
units of mass which constitute a nebula. The average square veloc- 
ity (21) might be derived from the shape of the spectral lines in the 
light from nebulae. Unfortunately, the practical determination of 
such shapes is at present exceedingly difficult, if not impossible. In 
addition the spectral lines in the light of nebulae are doubtless of 
complex origin, and the interpretation even of well-known shapes of 
lines is by no means an easy task. 


IV. NEBULAE AS GRAVITATIONAL LENSES 


As I have shown previously,° the probability of the overlapping of 
images of nebulae is considerable. The gravitational fields of a num- 
ber of ‘‘foreground”’ nebulae may therefore be expected to deflect the 
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light coming to us from certain background nebulae. The observa- 
tion of such gravitational lens effects promises to furnish us with the 
simplest and most accurate determination of nebular masses. No 
thorough search for these effects has as yet been undertaken. It 
would seem, perhaps, that if the masses of field nebulae were, on the 
average, as great as the masses of cluster nebulae obtained in sec- 
tion iii, gravitational lens effects among nebulae should have been 
long since discovered. Until many plates of rich nebular fields taken 
under excellent conditions of seeing have been carefully examined it 
would be dangerous, however, to draw any definite conclusions. 

The mathematical analysis of the formation of images of distant 
nebulae through the action of the gravitational fields of nearer nebu- 
lae will be given in detail in an article to be published in the Hel- 
vetica physica acta. 


V. STATISTICAL DISTRIBUTION IN SPACE OF DIFFERENT 
TYPES OF NEBULAE 

It will be shown elsewhere that the number of clusters of nebulae 
actually observed is far greater than the number that might be ex- 
pected for a random distribution of non-interacting objects. This 
tendency of nebulae toward clustering is no doubt due to the action 
of gravitational forces. 

By a bold extrapolation of well-known results of ordinary statisti- 
cal mechanics we adopt the following working hypothesis as a tenta- 
tive basis for the interpretation of future observations on the clus- 
tering of nebulae. 

BASIC PRINCIPLES 


1. The system of extragalactic nebulae throughout the known 
parts of the universe forms a statistically stationary system. 

2. Every constellation of nebulae is to be endowed with a proba- 
bility weight f(€) which is a function of the total energy ¢ of this con- 
stellation. Quantitatively the probability P of the occurrence of a 
certain configuration of nebulae is assumed to be of the type 


P= (£) (49) 
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Here V is the volume occupied by the configuration or cluster con- 
sidered, V, is the volume to be allotted, on the average, to any indi- 
vidual nebula in the known parts of the universe, and € is the total 
energy of the cluster in question, while €x will probably be found to 
be proportional to the average kinetic energy of individual nebulae. 
The function A(V/V.) can be determined a priori. On the other 
hand, f(€/€x) presumably will be found to be a monotonously de- 
creasing function in e/€x, analogous in type to a Boltzmann factor 


F = const X «v=. (50) 


Assuming the basic principles stated in the preceding to be cor- 
rect, we may draw the following hypothetical conclusions: 

a) The clustering of nebulae is favored by high values of f and is 
partially checked by low values of the a priori probability A. 

b) If, as would appear to be certain, nebulae are not all of the 
same mass, nebulae of high mass are favored in the process of clus- 
tering, since they contribute most to produce high values of the 
weight function f. 

c) As a consequence of b, we should expect that the frequency 
with which different types of nebulae occur will not be the same 
among field nebulae and among cluster nebulae. In other words, 
clustering is a process which tends to segregate certain types of 
nebulae from the remaining types. This may contribute toward the 
correct interpretation of the well-known fact that cluster nebulae 
are preponderantly of the globular and elliptical types, whereas field 
nebulae are mostly spirals. From the arguments put forth in the 
preceding section as well as in section iii it follows that it is not neces- 
sary as yet to call on evolutionary processes to explain why the repre- 
sentation of nebular types in clusters differs from that in the general 
field. Here, as in the interpretation of other astronomical phenom- 
ena, the idea of evolution may have been called upon prematurely. 
It cannot be overemphasized in this connection that systematic and 
irreversible evolutionary changes in the domain of astronomy have 
thus far in no case been definitely established. 

d) If cluster nebulae, on the average, are really more massive than 
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field nebulae, the conclusion suggests itself that globular nebulae 
may, somewhat unexpectedly, be among the most massive systems. 
It will be of great interest to check this inference by a search for 
gravitational lens effects among globular nebulae. 

The preceding considerations point toward the possibility of an 
entirely new approach in the study of masses of nebulae. We may 
argue somewhat as follows: 

The function A(V/V,), as said before, can be obtained from the 
theory of probabilities applied to random distributions in space of 
non-interacting objects. The function A, therefore, is known a pri- 
ori. The function f may be determined from counts of types of clus- 
ters of nebulae (single, double, triple, quadruple nebulae, etc.) in 
given volumes of space. Such counts will give directly the values of 
P characteristic for different clusters. The numerical values of the 
function f follow from f = P/A. In order to determine the masses 
of individual nebulae, we must express the argument ¢€/ex of the 
function f in terms of these masses and then seek to correlate each 
definite argument with one numerical value of f. To solve this prob- 
lem of the functional form of f we may proceed as follows: 

We first segregate the nebulae into classes of types 7;, 7.,73,.... 
T,. For these types the usual ones, E,, F,,.... E,, Sa, Ss, Se, etc., 
may, for instance, be chosen. As a first approximation we assume 
tentatively that the mass M,(L) of a nebula of a given type T, isa 
function of its luminosity L alone. The argument €/€x of the func- 
tion f may then be formulated mathematically in terms of M,, r,, 
and v,, where the M, are the unknown masses of the various types 
(T,) of nebulae in the cluster and where the velocities v, of these 
nebulae, as well as their mutual distances 7,,, are known from ob- 
servation. Since the numerical values of f are already known, the 
functional dependence of these values on the arguments €/€x can 
then be determined by a purely mathematical procedure. Once the 
form of the function f(€/€x) is known, the unknown masses follow 
from our knowledge of ¢/é€x expressed in terms of M,, 2,, 7,,. It 
should be noted, however, that the method just described gives only 
relative masses M,/M,, measured, for instance, with the mass M, 
of the type 7, taken as the arbitrary unit. Only if the absolute value 
of €x is known or if we have independent knowledge of M, can we 
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derive the absolute values M, from the statistics of nebular distribu- 
tion. 

Needless to say, this program, if it can be carried out, provides 
the most powerful method of determining the masses of all types of 
nebulae. In addition, it also enables us to determine the statistical 
weights f. The practical application of this method necessitates a 
great amount of observational work. In view of this fact, it will per- 
haps be advantageous to apply the preceding program first in a re- 
stricted form by a consideration of the distribution of various types 
of nebulae in one individual great cluster. The procedure to be ap- 
plied in this case is analogous to that used successfully by H. von 
Zeipel’ in his determination of the masses of different types of stars 
in certain clusters of stars. 

Since it is intended to carry out the investigation just mentioned 
on the Coma cluster, a few preliminary remarks concerning the dis- 
tribution of nebulae in this cluster are here given. 


VI. THE COMA CLUSTER OF NEBULAE 


According to Hubble and Humason,” the Coma cluster of nebulae 
“consists of about 800 nebulae scattered over an area roughly 1°7 in 
diameter. .... Photographic magnitudes range from about 14.1 
to 19.5, with 17.0 as the most frequent value.” The distance of 
the cluster is about 13.8 million parsecs. 

In order to get some preliminary data on the distribution of 
nebulae in the Coma cluster, photographs of this cluster were taken 
on Mount Palomar with the new 18-inch Schmidt-type telescope of 
the California Institute of Technology. The faintest nebulae which 
on limiting exposures (30-60 min) with this telescope can still be 
clearly distinguished from stars have an apparent magnitude close 
to 16.5. In Figure 3 dots represent nebulae which can be distin- 
guished on 30-minute exposures on panchromatic films. To avoid 
crowding of points near the center of the cluster, not all the nebulae 
in this region are marked which can be seen on the original photo- 
graphs. The counts given in Table 1 at different distances r from the 

10 EF. Hubble and M. L. Humason, Ap. J., 74, 131, 1931; see also the descriptions of 


the Coma cluster by M. Wolf, Heidelberg Pub., 1, 125, 1902, and H. Shapley, Harvard 
Bull., No. 896, 1934. 
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center of the Coma cluster include, however, all the nebulae which 
I have been able to identify on half-hour exposures. 

The nebula NGC 4874 (a 12"56™, 6 28°20’, 1930), which lies at 
(0°, o°) in Figure 3, was taken to be the approximate center of the 
cluster, no effort being made to determine a mathematically accurate 
central point. Concentric circles were then drawn around the adopted 
center, with radii r = mr,, where 7 is a whole number running from 
n = 1 to m = 32, andr, & 5 minutes of arc. The unit of area to 


TABLE 1 
COUNTS OF NEBULAE IN THE COMA CLUSTER 

63 1.799 1.20 .079 


which all counts are reduced is s = mré, or about 1/46.4 sq. deg. The 
numbers of nebulae per unit area, ”,, are averages for the ringlike 
areas which lie between r = nr, and r = (n + 1)ro. The first four 
figures in Table 1, however, are averages for the full circles the radii 
of which are r = 7o/5, 70/3, 7o/2, and ro, respectively. The corre- 
sponding numbers WN, of nebulae per square degree are NV, = 46.4 ,. 

In Figure 4 values of log, , are plotted against values of r. 
The full curve is drawn only approximately and does not correspond 
to any definite mathematical function. From the general character 
of this curve it is seen immediately that the Coma cluster extends to 
much greater distances than was originally assumed by Hubble and 
Humason.” At the edge of a circle the diameter of which is 4°5 in- 
stead of only 1°7, the average number of nebulae per unit area is 
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still higher than the corresponding number in the surrounding gen- 
eral field. Since our counts include only the brighter nebulae, it is 
to be expected that counts made with more powerful telescopes will 
enable us to follow the extensions of the Coma cluster still farther 
into the general field. 


log n, 


2.0°" T T T T T T 


o’ 25 50 75 100 125 150 175 7 


Fic. 4.—Counts of nebulae in the Coma cluster 


The high central condensation, the very gradual decrease of the 
number of nebulae per unit volume at great distances from the cen- 
ter of a cluster, and the great extension of this cluster become here 
apparent for the first time. It is quite as we should expect from the 
considerations of section v. According to these considerations, a 
cluster of nebulae analogous to an isothermal gravitational gas 
sphere may in some cases be expected to extend indefinitely far into 
space, until its extension is stopped through the formation of inde- 
pendent clusters in the regions surrounding it. 

The actual shape of the distribution curve in Figure 4 is also of 
great interest. We notice at once the great similarity of this curve 
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to the luminosity-curves of elliptical nebulae derived by Hubble.” 
According to him, the distribution of the intrinsic luminosity in 
globular nebulae corresponds very closely to the distribution of mass 
density in isothermal gas spheres as computed by R. Emden.’ The 
same is approximately true for the distribution of the brighter 
nebulae in the Coma cluster. This result also checks the general 
conclusions drawn from the basic principles which, according to the 
discussion in section v, determine the stationary configurations of 
clusters of nebulae. 

The total number ./ of nebulae in the Coma cluster that can be 
identified on photographs taken with the 18-inch Schmidt telescope 
is obtained as follows: The curve in Figure 4 apparently has a hori- 
zontal asymptote which corresponds to a value of m, not higher than 
N,, = 0.159. Thus, the corresponding number N of nebulae in the 
general field should be N = 46.4 X n,, = 7.38. The total number of 
nebulae counted to the distance r = 2°40’ from the center is 834. 
From this we must subtract 22.3 X 7.38 = 165 nebulae, since the 
area in question covers 22.3 sq. deg. The Coma cluster therefore 
comprises a number ./ of nebulae the brightnesses of which are 


greater than m = 16.5: 


Nw. = 670. (51) 


Since the most frequent apparent magnitude is about 17, we con- 
clude that .%,, is less than half the total number _/’ of nebulae in- 
corporated in the Coma cluster. This number must be at least equal 
to .V = 1500, and it may be even greater. 

Finally, a word with respect to the limiting magnitude m = 16.5 
which we have used in the preceding discussion. According to Hub- 
ble,? the number J, of nebulae per square degree which are brighter 
than m near the galactic pole is given by 


log Nm = 0.6m — 9.05. (52) 


Since the Coma cluster lies within a few degrees of the galactic pole, 
we may use relation (52) directly without introducing any correc- 
tion for partial obscuration. Our counts of nebulae lead to an aver- 


™ EF. Hubble, Ap. J., 71, 131, 1930. 
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age number of 7.38 nebulae in the general field surrounding the Coma 
cluster. Inserting N,, = 7.38 into the equation (52), we therefore 
find that m = 16.6 represents approximately the limiting magnitude 
at which it is still possible to distinguish images of average nebulae 
from those of stars on photographs taken with the 18-inch Schmidt 
telescope. 


VII. COMPARISON OF THE THREE METHODS 


Each of the three new methods for the determination of masses of 
nebulae which have been described makes use of a different funda- 
mental principle of physics. Thus, method iii is based on the virial 
theorem of classical mechanics; method iv takes advantage of the 
bending of light in gravitational fields; and method v is developed 
from considerations analogous to those which result in Boltzmann’s 
principle in ordinary statistical mechanics. Applied simultaneously, 
these three methods promise to supplement one another and to make 
possible the execution of exacting tests to the results obtained. 

Method iii can be applied advantageously only to clusters. Its 
application calls for the observation of radial velocities of cluster 
nebulae. The absolute dimensions of the cluster investigated also 
must be known. 

Method iv involves the observation of gravitational lens effects. 
Measurements of deflecting angles combined with data on the ab- 
solute distance of the “lens nebula” from the observer suffice to de- 
termine the mass of the lens nebula. The chances for the successful 
application of this method grow rapidly with the size of the available 
telescopes. Since method iii gives only the average masses of cluster 
nebulae and method v furnishes only the ratios between the masses 
of different types of nebulae, much depends on whether or not a sin- 
gle image of a nebula, modified through the gravitational field of 
another nebula, can be found. A single good case of this kind would, 
so to speak, provide us with the fixed point of Archimedes in our at- 
tempt to explore the physical characteristics of nebulae. 

Method v is the most powerful of all, since it enables us in princi- 
ple to find the masses of all types of nebulae, provided the absolute 
mass of a single type of nebula is known or that we have some inde- 
pendent way of finding a sufficiently accurate value of €x. Method v 
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also results automatically in the knowledge of the statistical weight 
functions f which govern the distribution of nebulae. The knowledge 
of these functions is of interest for two reasons: 

1. The weight functions derived from direct observations may be 
compared with those to be expected theoretically, for different ““mod- 
els” of the universe. Through such a comparison it should be possi- 
ble to decide whether the universe as a whole is in thermodynamic 
equilibrium™ *3 or is continually changing. 

2. It will be of particular interest, as proposed previously," to in- 
vestigate the probability P of the occurrence of clusters of nebulae 
in our “immediate” extragalactic neighborhood, as well as at great 
distances. If the universe is, for instance, expanding, we should ex- 
pect P to be different at different distances from the observer. The 
fact that a great cluster of nebulae, such as the Coma cluster, seems 
to represent a statistically stationary configuration suggests that a 
short time scale with 10° years as the characteristic age of the uni- 
verse is hardly adequate. Considerably longer time intervals would 
seem to be necessary to insure the formation of a stationary dis- 
tribution of nebulae in great clusters. A detailed analysis of the 
problem of the time scale, however, must be postponed until the dis- 
tribution of nebulae in a greater number of clusters has been investi- 
gated. 
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PHOTOELECTRIC MAGNITUDES AND COLORS OF 
EXTRAGALACTIC NEBULAE* 


JOEL STEBBINS' AND ALBERT E. WHITFORD” 


ABSTRACT 


The magnitudes of 165 extragalactic nebulae down to magnitude 15 have been 
measured with the photoelectric cell; of these, 112 have also been measured for color. 
The magnitudes, based upon the Harvard visual system, are reduced to the inter- 
national scale and give a correction of —o.1 to the photographic magnitudes of the 
Harvard Survey at magnitudes 11-13. 

The colors of the nebulae are referred to the color scale of dwarf stars. The nebulae 
of Hubble’s types E, Sa, and Sb are nearly uniform in color; those of type Sc are more 
varied, owing presumably to the presence of giant stars. 

There is a residual color excess for the nebulae, E= +-o0.14 mag. on the international 
scale, which is possibly due to selective absorption within the nebulae themselves. 
There is no dependence of color upon magnitude such as would be caused by selective 
absorption in intergalactic space. 


From the colors of nebulae, globular clusters, and B stars it is found that the ratio 
of selective to total space absorption is smaller in high than in low galactic latitude. 
The sun is assumed to be immersed in a dark layer of small selectivity near the galactic 
plane, while clouds of greater selectivity are present in nearly all directions along that 
plane but not immediately about the sun. 

The photoelectric cell furnishes a convenient means for comparing 
the total light of a diffuse object like a nebula with the light of a 
single star. The scale of distances for nebulae depends upon their 
assumed absolute magnitudes, which for near-by nebulae are found 
by comparing the relative magnitudes of the nebulae and the indi- 
vidual component stars. The photoelectric magnitudes of a dozen 
of the brightest nebulae have already been published by Whitford, 
who used the 10-inch telescope with focal diaphragms up to 30’ for 
some nebulae, but for the Andromeda nebula a 33-inch lens with a 
diaphragm more than 2° in diameter. The present work on fainter 
nebulae was begun by Stebbins in the summer of 1930. With a pho- 
tocell on the 60-inch reflector it was found possible to measure the 


* Contributions from the Mount Wilson Observatory, Carnegie Institution of Wash- 
ington, No. 577. 
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247 


248 JOEL STEBBINS AND ALBERT E. WHITFORD 


magnitudes and colors of the central areas of nebulae down to the 
eleventh or twelfth magnitude. The focal diaphragm was about 3’ 
in diameter for bright objects, and 1’ or less for faint ones. The work 
was continued with the same photometer on the 1oo-inch reflector 
in 1931 and 1932. From 1933 on, the amplifier was substituted for 
the electrometer in the installation. 

In addition to the nebulae, a few globular clusters were included 
in the program, and, after several clusters in low galactic latitude 
had been found to be strongly reddened, other clusters were added 
to the list; and, finally, all but two or three of the globular clusters 
observable from Mount Wilson were measured for color. The meth- 
ods of observation have been described in our paper on the clusters.‘ 
The measures and reductions of the clusters and nebulae were car- 
ried on together, so that the same system of magnitudes holds for 
both. The advantages of the amplifier over the electrometer were 
even greater for the nebulae than for the clusters, as the nebulae ex- 
tend to fainter magnitudes. 


COLOR SCALE 


With the same photocell and pair of filters there has been no 
change in the color scale except what was caused by variations of 
the mirrors of the large reflectors and by possible progressive changes 
in the cell over the years. For convenience the color scale was first 
determined from giant stars, but in the summer of 1936 about sixty 
dwarfs from the Mount Wilson Catalogue’ were observed to establish 
a scale comparable with the dwarf characteristics of the nebulae. 

The colors of the dwarf stars are in Table 1. The data in the first 
four columns are from the Draper Catalogue, except for the few faint 
stars marked C and 20C, which are from Cincinnati Publications, 
Nos. 18 and 20. The spectrum and the visual absolute magnitude M 
are from the Mount Wilson list. The last column gives the color in- 
dex C,, observed with our light filters and reduced to the standard 
season of 1932. We designate by C, the color index measured with 
the heavier filters used for the B stars.® Most of the stars in Table 1 


4 Mt. W. Contr., No. 547; Ap. - 84, 132, 1936. 
5’ Adams, Joy, Humason, and Brayton, Mt. W. Contr., No. 511; Ap. J., 81, 187, 1935. 
§ Stebbins and Huffer, Pub. Washburn Obs., 15, 218, 1934. 
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TABLE 1 
OBSERVED COLORS OF DWARF STARS 
HD R.A. 1900 Decl. 1900 Vis. m Spec. M C2 
+40° 23’ 8.7 Mo 8.9 +o™28 
© 29.4 —5 6 7.4 Go + .o1 
BOG +20 43 6.08 Ki 5-9 + 
Ce © 35.3 +39 39 7.52 K5 6.5 + .14 
© 39.9 + 1 15 8.14 Ks 6.5 + .16 
40.5 —13 25 6.11 F8 4-3 — .04 
PS, Ree © 43.1 + 4 46 5.82 K4 6.5 + .12 
(21s ee © 49.2 +86 47 8.9 G4 4.6 + .02 
2.8 + 1 28 6.69 G5 4.1 + .08 
I 14.0 — 9 27 8.9 Go 4-5 — .05 
+15 11 G6 4.0 + .I0 
+27 36 7.9 G6 4.6 + .03 
+42 7 5.10 Go 4.4 — .03 
TOOO?: 30:5 +19 35 6.23 G4 4.7 + .03 
I 42.9 +16 31 7.44 Go 2.8 + .o1 
2 9.2 +64 30 8.40 G5 4.9 .00 
2 25.4 +17 16 6.41 G5 + .13 
2 30.6 + 6 25 5.92 K4 6.7 + .16 
232.7 +30 24 7.21 Go 4.1 — .04 
2. 65.9 + 5 36 8.2 Ki + 
+ 8 27 K2 5-9 + .13 
72) 3 3200 +14 49 7.69 G4 4.2 + .05 
20000: +30 4.96 G5 4:7 
— § 42 8.1 K6 + .23 
$44.3 —20 9 8.2 Mo 8.4 + .27 
3:40:52 +22 23 7.8 G5 4.9 + .02 
14 30.2 +53 20 7 Ks 6.7 + .16 
+54 4 7.9 Ko 5-4 + .05 
OC 020.5555. 15 8.8 — 3 26 9.6 Mo 8.2 + .14 
TA75702 hes 16 16.5 +67 29 8.9 Mo 8.7 + .30 
+48 11 7.03 F8 — .06 
16: 33% —2 7 5.87 Ko 5-4 + .09 
16 41.4 +33 41 8.6 Mo 8.7 + .28 
16 47.9 + or 6.78 Go + .05 
16 53.8 +62 16 7.04 G4 + .or 
clr 2 16 59.8 +47 12 6.74 Ko 5.2 + .02 
Ch ee 16 59 8 — 454 7.90 Mo 7.8 + .20 
$60690. «06055. 18 6.3 +38 27 6.40 K2 5-9 + .09 
1006... 18 32.4 +45 39 9.8 M2 8.5 +o. 28 


| 
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TABLE 1—Continued 
HD R.A. 1900 Decl. 1900 Vis. m Spec. M C2 

19"29™5 +58° 23’ 6.70 K5 6.1 +o™10 
TOSIOZ... ss 20 13.8 +76 55 9.3 Mo 8.9 + .31 
ar + 6 41 8.9 K6 + .23 
2, 21 8.8 +73 18 8.8 K4 6.8 + .06 
21 44.1 + 5 15 8.61 K6 + .16 
22 4.2 — 8 2 6.63 Go + .03 
22 30.5 + 4 52 8.4 G4 4.4 .00 
—o 6 8.0 G4 4.0 + .09 
| rr 22 46.4 +13 26 8.0 K4 6.4 + .08 
SS ae 22 47.3 +31 12 9.4 K6 6.9 + .17 
22 52.6 +20 14 5.59 Go 4.3 .00 
— 4 23 7.60 K4 6.1 + .08 
23 17.8 — {I 8.1 Kr + .I0 
23 26.4 — 4 38 6.50 F8 3.6 — .07 
ei Os 23 30.4 +30 27 6.72 Go 4.4 — .04 
23-4120 +17 53 8.0 G5 4.9 + .o1 
232500. 23.-30..2 +010 8.03 F8 453 — .05 
23 40.0 +29 0 8.9 K2 + .10 
23 44.0 + 1 52 9.1 M2 9.2 +0. 25 

C; 

—0.3 

— 

I 
oF 

+ 

+ .2- 

+ .3-L 

A F G K M 
Fic. 1.—Colors of dwarf stars. Heavy line, mean for dwarfs; broken line, mean for 


i 
| 
giants. | 
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were observed only once. Because of the dispersion of color within a 
spectral class it seemed better to take additional stars than to repeat 
the measures of each star. The colors C, in Table 1 are shown in 
Figure 1, and the smoothed means are in Table 2. There is some 


TABLE 2 
NORMAL COLORS OF GIANT AND DWARF STARS 
Cotor InpDEx, C2 Cotor Inpex, C2 
SPECTRUM SPECTRUM 

Giants Dwarfs Giants Dwarfs 
—0.02 


evidence of an absolute-magnitude effect within the spectral classes, 
particularly in G4, G5, and G6. In Gs the dependence of C, upon M 
is pronounced: 


M M Cs 


In other spectral classes the data are insufficient to show such a rela- 
tion; but a rough attempt was made to find the color for the mean 
absolute magnitude of each spectral subdivision, and the curve in 
Figure 1 was drawn with some allowance for this effect. In Table 2 
the previously derived colors for giants are also included. 

We are well aware of the compressed scale of our colors, which has 
a difference of only 0.34 mag. between Ao and Ko for dwarf stars. 
We found from the Polar Sequence that the relation to the interna- 
tional scale is given by 


IC = 2.62 C, + Const. 
For bright stars we get a probable error of less than + 0.01 mag. for 


one observation of color, corresponding to + 0.02 mag. or +0.03 mag. 
on the international system. The errors for faint nebulae are, of 
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course, larger. The differences in color between stars with the same 
spectral lines are already conspicuous on our scale; a combination of 
cell and filters using more widely separated spectral regions would 
simply exaggerate these differences in the stars themselves. 


STANDARD MAGNITUDES 
The photoelectric magnitudes of the nebulae were referred to the 
Harvard visual magnitudes of stars, mostly between the sixth and 
the seventh magnitude. The observations of a night usually began 


TABLE 3 
MAGNITUDES OF STANDARD STARS 
Group HD Spec. HV C; Pe Pey |Pe—Pey| Resid. 
(1) (2) (3) (4) (5) (6) (7) (8) (9) 
149631} 7.30 | —oM23] 7.30 7.37. | —oMo7 | —oMo6 
149632} Ao 6.27 — .24 |] 6.43 6.32 + .11 | + .12 
Fa2 6.32 — .12] 6.63 6.57 | + .06 | + .07 
156717 Ao 6.04 — .23 6.01 6.11 = ito | — 360 
160018} Ko 5.92 | +0.26 |] 6.75 6.81 —0.06 | —0.05 
Group Pe — Peg Resid 
—oMor —oMo4 
.00 — .03 
+ .or — .02 
+o.11 +0.08 
Mean..... +0.03 +0.04 


with several standard stars measured in the twilight, and other 
standards were taken from time to time between the measures of 
the nebulae. Usually the conditions were so good that all the stand- 
ards observed over several hours could be used as a group; but when 
there were interruptions, the measures and reductions of the night 
were broken into several sections. The reductions of the standard 
stars from visual to photoelectric magnitudes were determined by 
the measured colors. Since the 1oo-inch cannot be pointed near the 
pole, a comparison of some of the standard stars with the North 
Polar Sequence was made with the 60-inch reflector. As it is not con- 
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venient to rotate the Newtonian cage of the 60-inch during the night, 
the southern stars were compared with the Polar Sequence by way of 
intermediate stars culminating near the zenith, which could be 
reached with the cage in either the north or the south position. In 
the tests were two southern groups, I and II, of five stars each, and 
two intermediate groups, III and IV, of five and six stars, respec- 
tively. 

In Table 3, HV is the Harvard visual magnitude, C, the photo- 
electric color, Pe the observed photoelectric magnitude referred to 
our arbitrary zero point of the Polar Sequence, and Peg the photo- 
electric magnitude computed from the visual magnitude by the 


formula 
Pey = HV + 1.68(C, + 0.27). (1) 


The mean of Pe — Pey for each group is the difference between the 
zero points from the Polar Sequence and the Harvard visual mag- 
nitudes. Taking the average for the four groups, we have 


Pe = Pey + 0.03. (2) 


For the polar stars we found 


Pg, = Pe — 0.05 + 0.53(C2 + 0.21) , (3) 


where Pg, is the photographic magnitude computed from the photo- 
electric. The constant —o.o5 brings Pg, into agreement with J Pg 
for the mean of five stars at magnitude 6.0. Substituting (2) in (3), 
we have 


Pg, = Pew — 0.02 + 0.53(C2 + 0.21). (4) 


The constant term —o.02 has been determined here from only 21 
stars, and it must vary with any systematic errors of the Harvard 
visual magnitudes over different parts of the sky. The term can be 
discarded, and we adopt 


Pgy = Pe+o.53(C. + 0.21), (5) 


where Pe now has the same zero point as Peg for any group of stars. 
This close agreement in zero point between the Harvard visual 
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scale and the international photographic scale is accidental. In equa- 
tion (1) we adopted C, = —o0.27 for an Ao star in a nonobscured re- 
gion; but because of the selective absorption at the pole, C, = 
—o.21 when the international color index JC = 0.00. To keep the 
magnitudes on the international system throughout, we should use 
0.21 in place of 0.27 in (1), which would change the constant term in 
(4) by 0.10. Instead, we let the photoelectric magnitudes remain as 
they are, and then the photographic magnitudes derived from (5) are 
nearly enough on the international system. On our scale the zero 
point of Pg, is made to agree with Pg at magnitude 6.0. The preci- 
sion of the magnitudes of the comparison stars is indicated by the re- 
siduals in column (9) of Table 3. Usually the average residual for 
a star islessthan +o0.10 mag.; the means from groups of several stars 
are good enough for standards for the nebulae. 

It might appear better to refer the magnitudes of nebulae to 
standards in Selected Areas than to the Harvard visual magnitudes 
of bright stars, but the Harvard stars are easily picked up and are 
available all over the sky. The manipulation of so large a telescope 
as the 100-inch is a slow process, and we could not afford to spend 
valuable time in identifying and checking faint standard stars. From 
laboratory and other tests in observing the North Polar Sequence 
with the 60-inch, we found the scale of magnitudes with the photo- 
cell to be correct within o.o1 mag. or 0.02 mag. over the range from 
the second to the twelfth magnitude. Hence the same scale would 
hold to the thirteenti: magnitude with the 1oo-inch, and the exten- 
sion to one or two magnitudes still fainter is a small extrapolation. 


THE OBSERVATIONS 


The present results for extragalactic nebulae are brought together 
in Table 4. Column (1) gives the New General Catalogue number; 
column (2), the year of observation (measures after the aluminizing 
of the large mirrors in 1935 are marked “35a’’). The photoelectric 
magnitudes in columns (3) to (6) are for the respective diaphragms. 
Usually the 64” diaphragm was used first for both the magnitude 
and the color; then the measures through other diaphragms were 
made differentially from 64”. In column (7) all measures on the 
same night were counted as one observation. In column (8) the com- 
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puted photographic magnitude Pg, was derived from the Pe mag- 
nitude for the largest diaphragm by equation (5). In cases where 
the color C, had not been observed, the mean C, for the type of 
nebula was used. 

In column (9) is the Harvard photographic magnitude by Shapley 
and Ames.’? Where our diaphragm was clearly large enough to in- 
clude all of the nebula indicated by the recorded diameter, the dif- 
ference Pg, — HPg is entered in column (10). The diameters in 
column (11) were taken from the Harvard Survey, with some addi- 
tions or changes furnished by Hubble from Mount Wilson plates. 
The types of the nebulae in column (12), as well as the spectra in 
(13), were also furnished by Hubble from plates usually taken by 
Humason. 

In column (14) the measured color index C, refers to the dia- 
phragm 64”, or to the smallest of the other diaphragms if 64” was 
not used. The seasonal corrections which have already been applied 
to C, in Table 4 are: 1932, 0.00 mag.; 1933, —0.01 mag.; 1934, —0.05 
mag.; 1935, —0.05 mag.; 1935a, to.05 mag. The jump between 1935 
and 1935a was due to the aluminizing of the mirrors of the telescopes. 

The average deviation and the number of nights are in columns 
(15) and (16), while in (17) and (18) are the approximate galactic co- 
ordinates taken from the Harvard Survey or the Lund Tables. The 
galactic pole is at R.A. 12"40™, Decl. +28° (1900). 

In the summer of 1936 we measured a dozen nebulae in low lati- 
tudes near the edge of the zone of avoidance. These nebulae are 
faint, of magnitude 14 or 15, and are likely to be involved in rich 
star fields; hence small diaphragms were used to avoid including 
stars and to reduce the light of the sky. Measures were not attempt- 
ed if there were superimposed stars. The rate-of-drift method of 
measurement was used, with photographic registration of the gal- 
vanometer record. The recording drum was driven at the rate of 
about 50 mm per minute, and the slope of the traced line was deter- 
mined for an interval of 30 or 45 seconds. The results are in Table 5. 

Magnitudes.—The magnitudes for the different areas of the nebu- 
lae in Table 4 indicate both the advantages and the disadvantages 
of measures with the photocell. The limited diaphragms tell exactly 


7 Harvard Ann., 88, No. 2, 1932. 
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which part of a nebula was measured, but for large nebulae the mag- 
nitude of a central area may not be of much significance. For in- 
stance, for a diaphragm of 42” the photoelectric magnitudes of the 
Andromeda nebula and its brighter companion, NGC 224 and 221, 
are 9.33 and 10.03, respectively. The total magnitude of the com- 
panion is reached at about 9.33 for 163”, but it is necessary to open 
to a diameter of several degrees to get most of the main nebula, 
which, with allowance for the many superimposed stars, turns out 
to be of magnitude about 4.5. It was not feasible to make a special 


TABLE 5 
MEASURES OF FAINT NEBULAE, 1936 

NGC Pe Pg» | Diaph.| Diameters | Type | Spec. C2 I b pa 
13.41| 13.5| 32” | 80” X4o"|Ir |Gpec.| +oMos5} 119°|—12° 
¥4.74| 14.9) 32 |'25 X.. IB3 [G3 |—12 
TOBE TA: 00} 1552) 32 .I9] 119 |—12 
38.7) 15 | .18| 11 |-+12 | 2 obs. 
14.36] 14.5] 24 | 90 X20 |Sa_ |G3 .14] 19 |-+12 
14.01] 14.2} 24 |95 X65 |Sa_ |G5 |-+12 


study of every nebula on the list and be sure that we had all the 
outer regions. The magnitudes, as given, furnish a series of stand- 
ards which can be used for a further study of photographs to get 
over-all magnitudes if desired. 

While the magnitudes depending upon the comparison stars may 
have probable errors up to about +0.1 mag., the values for the four 
diaphragms in Table 4 are given to hundredths because their dif- 
ferences are significant to two decimals. There are numerous cases 
where the difference between, say, 128” and 163” indicates faint 
outer extensions which are not always shown on photographs. The 
cases where these extensions are real seem to run from about mag- 
nitude 24 to 25 per square second of arc. The light of the dark non- 
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galactic sky at Mount Wilson is something like magnitude 22 per 
square second. Most of the measures with different diaphragms on 
the same nebula were made in the early years with the electrometer; 
the later measures with the amplifier were on fainter and smaller ob- 
jects. Under the most favorable conditions the limit of detection of 
faint nebulosity is from magnitude 26 to 27 per square second, or 
about 4-5 mag. fainter than the surface brightness of the sky. Such 
a limit would be independent of the size of the telescope were it not 
for the interference of the stars. However, with a large telescope the 
number of faint stars brought into view is not sufficient to compen- 
sate for the reduction of the field, and it is therefore easier to find a 
blank field with the too-inch than with a smaller telescope of the 
same focal ratio. 

The photographic magnitudes Pg, in column (8) of Table 4 are 
directly comparable with the Harvard magnitudes in (9) when the 
diaphragm used was as large as the major diameter in (11). In these 
cases, where our measure took in all the nebula, the difference 
Pg, — HPg is given in column (10). The mean value of the differ- 
ence Pg, — H Pg for 59 nebulae is —0.08 +0.02(P.E.),which happens 
to agree exactly with the difference Yerkes minus Harvard found 
by Keenan* from photographic magnitudes of about fifty nebulae 
between the eleventh and thirteenth magnitude and north of declina- 
tion +50°. Seyfert? had also found —o.05 for the correction to the 
Harvard magnitudes fainter than 11.0. Since our magnitudes were 
connected with the Polar Sequence by way of the Harvard visual 
system, the close agreement of the three results must be accidental. 
From our determination Hubble has applied a correction of —o.1 to 
the Harvard magnitudes at magnitude 13. 

The mean error of a single difference Pg, — H Pg is + 0.25, whichis 
about the same as the + 0.26 found by Keenan for Yerkes minus Har- 
vard. The internal accidental error of Pg, is small compared with 
the errors of the Harvard visual magnitudes, and it seems fair to 
assume a mean error of +0.15 for Pg, and +0.20 for HPg. Our 
measures show that the Harvard magnitudes of nebulae from 11 to 
13 are better than the ordinary Durchmusterung magnitudes of 

8 Ap. J., 85, 326, 1937. 

9 Harvard Circ., No. 403, p. 10, 1935. 
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stars, and that they make a reliable basis for magnitudes of fainter 
nebulae. . 
Colors.—The color indices in Tables 3 and 4 are independent of 
the magnitudes. With C, referring to our heavy filters and C, to the 
light filters used for the nebulae, the relation between C, and C, is 


C, + 0.10 = 1.68(C, + 0.27), (6) 


where the zero point for each system is fixed by the average color 
index of an unobscured Ao star. The relation between the color 
excesses in the two systems is, therefore, 


E, = 1.688, . (7) 


The colors of some nebulae were not measured because they were 
too faint and difficult or because they would have taken more time 
than we could afford to spend on the observing program. The mag- 
nitude of a nebula from a couple of galvanometer deflections is 
easily determined as accurately as the magnitudes of the comparison 
stars; but the colors must be more precise, and hence are more diffi- 
cult. The errors of the colors increase, of course, with the faintness 
of the nebulae. From observations of the same nebulae on different 
nights the mean error of one observation of C, for different mag- 
nitudes is as follows: 


Brighter than magnitude 12.0............. +0044 
Fainter than magnitude 12.0.............. .070 


The light of even the darkest sky adds up to about magnitude 13 for 
a 64” circle; and when the superimposed sky is as bright as the nebu- 
la, the observations are difficult by any method. With the amplifier 
and with alternate readings on nebulae and blank sky through small 
diaphragms, colors can be measured down to magnitude 15, as in 
Table 5; but the observations are tedious and require about an hour 
for each object. 

The study of the colors of extragalactic nebulae was originally 
undertaken to determine how much selective absorption there might 
be either within our own galaxy or outside in intergalactic space. 
Perhaps the first test is to compare, for each type of nebula, the ob- 
served mean color with the corresponding color for the mean spec- 
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tral class, as in Table 6. The mean spectral classes for the different 
types of nebulae are given by Hubble.” For each type the mean col- 
or index C, for dwarf stars of the same spectrum is taken from our 
Table 2, and the difference E, is the mean color excess. The mean 


TABLE 6 
MEAN COLORS OF NEBULAE 
Spectral Mean Cz |Observed C2 M.E. 

Tepe Class Stars Nebulae 1 Neb. 
G4 +o™o2 | +oMog | +0™Mo7 | +0Mo53 30 
| ees G3 + .o1 + .07| + .06 .046 13 
G2 .00 | + .06| + .06 .053 12 
Fo —0.04 | —0.06 | —0.02 | +0.078 36 


error for a nebula is computed from the observed deviations of C, 
from the mean for each type. For this test, nebulae in latitudes less 
than 20° were excluded. 

One defect of the material in Table 6 is that the mean spectral 
classes and the observed colors refer only in part to the same nebulae. 


TABLE 7 


MEAN COLOR EXCESSES OF NEBULAE 
WITH KNOWN SPECTRA 


M.E. 
1 Neb. No 
+006 +oMos1 13 
+ .07 .043 8 
+ .06 .058 7 
—0.04 +0.099 7 


In Table 7 are the mean color excesses of those nebulae for which 
we have both the individual spectra and the colors. These data are 
included in Table 6; but on any basis of computation the three 
types E, Sa, and Sb show a color excess E, of +0.06 mag. or +0.07 
mag., with a probable error of +o0.o1 mag. to +0.02 mag. for each 
mean. The Sc nebulae have earlier spectra and are bluer in color. 
Seares™ found in a number of late-type spirals that the outer arms 
are distinctly bluer than the central regions, but even in the central 


%© The Realm of the Nebulae, p. 53, 1936. 
= Mt. W. Comm., No. 36; Proc. Nat. Acad. Sci., 2, 553, 1916. 
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regions the presence of a few white or blue giant stars would cause 
these nebulae to appear bluer than the other types of nebulae whose 
spectra seem to be predominantly dwarflike. Presumably, there 
could not be a large proportion of white or blue giant stars in a nebu- 
la without their showing up in the spectrum, but we must have some 
leeway in interpreting the integrated spectrum and color of a million 
stars. 

The mean error of one observed C, from the measures of the same 
nebulae on different nights was +0.053 mag., which, compared with 


TABLE 8 

COLOR EXCESS AT DIFFERENT LATITUDES 
b csc 5 E; (Obs.) | Ex (Comp.) Oo-C No. 
(1) (2) (3) (4) (5) (6) 
9°—17° 4-7 +o™17 +o™17 oMoo II 
20 —27 2.6 53 .13 .00 10 
Nebulae. .... 32 —49 .10 II — .o1 15 
54 —69 .08 — .03 12 
73 —76 .10 + .03 18 
20 —23 2.8 05 04 + .o1 3 
Globular clus-| | 28 —34 2.0 02 03 5 
37 —50 .03 .03 00 6 
73 —89 1.0 .02 .02 .00 5 
8 —14 .O4 .00 10 
15 —19 .02 03 — 15 
20 —27 2.6 .02 .02 .00 18 
30 —37 1.8 03 .02 + .o1 14 
40 —62 +0.00 +0.o1 13 


the mean errors or standard deviations in the next to the last columns 
of Tables 6 and 7, shows that for the types E, Sa, and Sb the 
dispersion in color is practically accounted for by errors of observa- 
tion. The greater dispersion for the Sc nebulae is probably a real 
effect, caused by varying proportions of giant stars. 

The colors of all the observed nebulae of types E, Sa, and Sb are 
grouped according to galactic latitude in Table 8.% The respec- 
tive columns contain: (1) the limiting latitudes; (2) the mean csc b 


2 The single nebula with a large color excess, IC 356, E2= +0.39 mag. for b=14), is 
not included in the mean. This nebula was difficult with the 60-inch, but the reddening 
seems to be real. The near-by nebula, IC 342, in lower latitude has E,=-+0.12 mag. 

The globular clusters NGC 2419 and 7492 were also omitted as being anomalous; 
and NGC 6779, as being on the border of the zone. 
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for the group; (3) the mean color excess on the scale of E, for com- 
parison with globular clusters and B stars, E, = 1.68E,; (4) the com- 
puted E, from equations (8); column (5) gives the residual, and (6) 
the number of objects. Since the observed E, is derived by sub- 
tracting the color for the mean spectral class from the measured color 
of each nebula, the difference is not necessarily zero for high latitude. 


GLOBULAR CLUSTERS AND B STARS 


The nebulae can be compared with the globular clusters and B 
stars outside the nebular zone of avoidance. The clusters are prob- 
ably all far enough from the galactic plane to be well outside the 
hypothetical absorbing layer, while the B stars are those with com- 
puted distances from the plane, 7 sin b, > 100 parsecs. The mean 
color excesses of these clusters and B stars are also given in Table 
8; and the solutions for their variation according to the cosecant 
law give equations (8), with probable errors. 


Nebulae, E, = +0“o88 + o™or7 csc b 


x 
Clusters, E, = +0.009 + 0.013 csc b (8) 
+16 


B stars, EF, = +0.004 + 0.007 csc b 
+ 6 + 4 


In these equations the three constant terms give the computed color 
excess for zero absorption. For the nebulae the term may be real. 
For the clusters and B stars E, should equal 0.00 mag. for csc b = 1; 
the small discrepancies are due to the choice of the zero points. The 
three coefficients of csc b indicate a small but probably real variation 
with the latitude in each case. Some of the B stars may still be in 
the absorbing layer, and hence show a smaller latitude effect. 

We can adopt the mean from the nebulae and clusters as the selec- 
tive absorption at the pole, EH, = +0.015 mag. +0.005 mag. (P.E.), 
which is in contrast with Hubble’s value of 0.25 mag. for the total 
photographic absorption. 

These figures indicate that in high latitudes the selective absorp- 
tion is a small fraction of the total. If we assume the same propor- 
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tion of selective absorption in the zone of avoidance, we run into 
difficulties. There are five strongly colored globular clusters which 
average E, = +0.72 mag.; the total absorption would then be 0.72 
mag. X 0.25/(0.015+0.005), or 12 mag., with a probable error of 
+4 mag. Similarly, the many B stars with values of EZ, between 
+o.4 mag. and +0.6 mag. would require an impossibly large total 
absorption. Despite the uncertainty of the numerical results, there 
seems to be no way of reconciling the small selective absorption in 
high latitude and the large selective absorption in low latitude with- 
out abandoning the assumption of a constant ratio of the selective 
to the total absorption. In the course of various discussions with 
Dr. Hubble he has suggested that the space absorption in the Galaxy 
may involve a tenuous layer of small selectivity, condensed clouds 
of large selectivity, and the source of stationary spectral lines. Our 
material can be made to agree with this view by assuming that the 
sun is in a region of small selectivity. The absorbing stratum ob- 
scures the nebulae roughly according to the cosecant law; and 
when the line of sight approaches the galactic circle, it passes 
through the clouds of large selectivity, but the nebulae are then lost 
to view. The nebulae, clusters, and B stars all agree in showing that 
both the selective and the total absorption are spotted in latitude 
and longitude. The cosecant law can be only a rough approxima- 
tion, as is shown by the irregular outlines of the zone of avoidance. 
Even the naked-eye appearance of the Milky Way contradicts a 
simple geometrical conception of the Galaxy. 


INTRINSIC COLORS OF NEBULAE 


If the selective absorption within our Galaxy has little effect upon 
the color of the visible nebulae, we can study the intrinsic colors of 
the nebulae themselves. The constant term in the first equation (8), 
E, = +0.088 +0.007 mag., is the unexplained color excess in the 
nebulae. The natural interpretation of such a result is that it can 
be due to some systematic error. In terms of spectral class and 
color the discrepancy means that nebulae of estimated spectrum G3 
have the color gg. If the nebulae of types E, Sa, and Sb contained 
giants and dwarfs intermixed, instead of dwarfs only, as assumed, 
the difference would be less but would not vanish. However, the 
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giant stars are not found in near-by nebulae like M32, whose spec- 
trum is certainly dwarflike. 

The color excess does not seem to be a function of the nebular dis- 
tance, but we have only a small range of magnitude for the test. In 
Table 9 the 30 type-E nebulae are divided into three groups, ac- 
cording to apparent magnitude. There is no indication of increasing 
color index with increasing distance; the small reverse effect must be 
accidental. The nebulae of magnitude 12.7 are probably ten times 
as far away as M32, but their average color is about the same. This 
test of the transparency of space can be continued outward until the 
red shift begins to affect the colors of the nebulae. 


TABLE 9 
m Ex P.E. No. 
II .024 10 
+0.08 +0.017 10 


We are left with the residual color excess, E, = +0.09 mag., as 
a possible measure of selective absorption within the nebulae them- 
selves. It is simple to assume an arrangement of scattering material 
within a nebula sufficient to produce this effect for an external ob- 
server. A thin layer of dark stuff near the main plane would redden 
all the stars on the far side of a nebula. Our own Galaxy may well 
have something of the sort near the nucleus, but near the sun the 
effect at right angles to the main plane would be small. Because of 
the manner of finding the color excess of the nebulae, we attach little 
weight to the numerical value; but the discrepancy between esti- 
mated spectrum and measured color is in the right direction, to be 
accounted for by selective absorption within the nebulae. 


SUMMARY 


We summarize the main conclusions from the colors of B stars, 
globular clusters, and nebulae. In addition to the published re- 
sults, * ° we have the colors of more than four hundred B stars, in- 
cluding most of those in classes Bo to B2 north of declination — 40° 
and down to the limit of the Draper Catalogue. For uniformity the 
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colors are all expressed in terms of E,; the relation to the internation- 
al scale is JE = 1.55E;. 

The brightest, and hence the nearest, B stars are all about normal 
in color, Z, = o for all latitudes. The fainter B stars in the flattened 
system are near the galactic circle. They show increasing color index 
with increasing distance from the sun, and practically all stars with 
E, > +0.10 are located within the zone of avoidance of the nebulae. 
Outside this zone the B stars show a small latitude effect, represented 
roughly by E, = (+0.007 +0.002) csc }, from stars more than 100 
parsecs from the galactic plane. The reddest B stars are usually 
within 3° or 4° of the galactic circle; they have color excesses up to 
E, = +0.30 in all longitudes from 300° through 360° to 200°. The 
reddening is rather spotted along the Milky Way, but in general it 
is greatest near the galactic center and least toward the anticenter. 
Color excesses of E, = +0.63 and +0.68 have been found for two 
eighth-magnitude Bo stars in longitude 345°, with a dozen others 
near by with EZ, = +0.40 to +0.55, but we have not found a single 
B star with EZ, > +0.30 in the range of longitude from 115° to 190°. 
There are, however, not many B stars of magnitude 8.0 and fainter 
in this region of the anticenter. The coefficient of selective absorp- 
tion—color excess per kiloparsec—is perhaps twice as great toward 
the center of the Galaxy as toward the anticenter. The reddening of 
the B stars is not primarily connected with obviously obscured 
areas; in fact, the brightest star clouds in Sagittarius contain many 
reddened B stars. Their colors are likewise not directly connected 
with the interstellar absorption of calcium and sodium. However, 
with Dr. Merrill, we have a suspicion of a closer connection between 
our color excess and the strength of the new unknown interstellar 
lines he has found in the spectra of early-type stars. 

The colors of the globular clusters are as strikingly related to the 
galactic center as are the positions of the clusters themselves. A cir- 
cle with about 35° radius and its center at longitude 330° includes 
some 54 globular clusters, more than half of those known. Within 
this circle, which also contains the great bulge of the zone of avoid- 
ance of the nebulae, are all the clusters with marked reddening. The 
5 reddest clusters, E, = +0.64 to +0.82, are on the edges of the 
middle zone, 7° wide, where no globular clusters can be seen. The 
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22 clusters outside the central group, which we have measured, are 
widely scattered over the sky and are all nearly of the same color. 
The small latitude effect is represented by EZ, = (+0.013 + 0.009) X 
csc 6. The boundaries of the nebular zone of avoidance are not sharp. 
The globular clusters are the redder the fewer the nebulae in the 
field, until, when E, > +0.20 for a cluster, the nebulae are blotted 
out. This limit seems to be about the same for B stars. 

The nebulae of types E, Sa, and Sb are nearly uniform in color; 
those of type Sc are more varied, owing presumably to the presence 
of giant stars. The colors of the first three types indicate an excess 
of E, = +0.09, possibly because of selective absorption within the 
nebulae. There is no indication of selective absorption in inter- 
nebular space. The latitude effect for the nebulae is represented by 
E, = (+0.017+0.005) csc b. Thus the nebulae, like the B stars and 
the globular clusters outside the zone of avoidance, show little varia- 
tion in color. Despite the paradox, what we need are colors of nebu- 
lae within that zone. There are a few “‘holes”’ in low latitude where 
very faint nebulae are detected, and we hope to observe some of 
these if any as bright as the fifteenth or sixteenth magnitude can be 
found. 

In order to reconcile the small selective absorption in high latitude 
with the large selective absorption in low latitude, we assume that 
the sun is immersed in a layer of small selectivity near the plane of 
the Galaxy. This layer would account for the total photographic 
absorption shown by the nebulae, which varies roughly according to 
the cosecant law and amounts to 0.25 mag. at the pole. The clouds 
of large selectivity which produce the strong reddening in low appar- 
ent latitude are also near the plane of the Galaxy but not immediate- 
ly about the sun. They cannot be far distant along the main plane 
since some of the B stars are strongly colored at only a few hundred 
parsecs. 

In the study of the nebulae the color effects we were looking for 
have all turned out to be very small; but the side results of conspicu- 
ous reddening of B stars and globular clusters toward the galactic 
center, and the lesser reddening of B stars along the galactic circle 
in all longitudes, have indicated where and how much the nebulae 
may contribute to the problem. Because of the uniform intrinsic 
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colors of the nebulae, particularly of type E, it is hoped that an ex- 
tension to fainter objects will produce results barely hinted at in the 
present investigation. 


We are indebted to Dr. Edwin Hubble for constant collaboration 
and advice. This investigation has been supported by grants from 
the Alumni Research Fund of the University of Wisconsin, the Na- 
tional Research Council, and the California Institute of Technology, 
the last in connection with plans for the 200-inch reflector. 
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ABSTRACT 


The main purpose of this Contribution is to record measurements of the intensities 
and displacements of the interstellar lines D1, D2, H, and K in the spectra of numerous 
early-type stars. Intensities of the unidentified lines \X 5780, 5797, and 6284 are in- 
cluded. Observational methods and the technique of measurement are fully described, 
and the accuracy of the results is discussed. 

Comparison of the Mount Wilson values of the intensities of the D lines with those of 
Beals (Table 3, Fig. 1) indicates a fairly small accidental error in both series but a con- 
siderable systematic difference between them. Displacements of the D lines measured 
at various observatories (Table 7) are in good agreement. 

The Catalogue, Table 6, includes 254 stars brighter than apparent magnitude 7.0, 
147 fainter stars, and, in addition, 3 novae. The numbers of objects with the various 
kinds of data are as follows: 


Ca Na dA 6284 
Displacement... .. . 174 264 
: measured........... 114 213 129 
Intensity estimated.......... 89 107 


Table 8 lists stars with especially strong interstellar lines. 


PROGRAM 


When observations of interstellar lines were begun at Mount 
Wilson many years ago, the plan was to concentrate on a few small 
regions of the sky in which spectra of all B-type stars down to the 
lowest practicable limit of magnitude were to be examined. It was 
hoped in this way to determine the behavior of detached lines as a 
function of distance from the sun. 

As the general complexity of the whole problem of the intensities 
and displacements of these lines became increasingly evident, how- 
ever, the program was enlarged to include numerous stars in other 
regions. One important consideration was the desire to provide ade- 
quate material for a systematic investigation of interstellar lines in 
relationship to galactic rotation. This necessitated including stars 


* Contributions from the Mount Wilson Observatory, Carnegie Institution of Wash- 
ington, No. 576. 
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near the galactic plane at all longitudes within reach from Mount 
Wilson. Another extension of the program arose from the emphasis 
placed, in the latter part of the work, on stars of high luminosity. 
Stars of types earlier than B2, and especially c stars of types Bo to 
A2, offer substantial strategic advantages in extending the observa- 
tional limits outward into space. A c star (type Bo to A2) of the 
seventh apparent magnitude is probably as distant as one of type 
B8n of the twelfth magnitude. 

The distribution in longitude and latitude of the stars finally ob- 
served resembles roughly that of B-type stars in general, except for 
decided concentration in two small areas, one in Cygnus, / = 40°, 
b = +1°, the other surrounding the h and x clusters in Perseus, 
1 = 102°, b = —3°, and less conspicuous concentration in a larger 
area in Cepheus, / = 70°, b = +5°, and another in Cassiopeia, / = 
85°, = o°. 


OBSERVATIONS OF LINES IN THE YELLOW AND RED 


Spectrograms for the D lines and A 5780, 5797, and 6284, taken 
with a plane-grating spectrograph" at the Cassegrain focus of the 
100-inch telescope, number approximately 720. Of these, about 660 
were taken between May, 1928, and August, 1936, with the instru- 
ment in its specially designed mounting; the remainder were taken 
between September, 1924, and March, 1926, with the grating in a 
preliminary mounting. 

The majority of the spectrograms were taken with an 18-inch 
camera, dispersion 34 A/mm; the remainder with a 10-inch camera, 
66 A/mm. Most of the spectra were widened. Comparison lines of 
neon have been employed throughout, while the D3 line of helium 
has been included since October, 1932, and the yellow mercury lines 
since May, 1936. 

Since 1930, the III C emulsion supplied by the Eastman Kodak 
Company has been regularly used. This fine-grained, contrasty emul- 
sion is very rapid when hypersensitized by bathing in dilute am- 
monia, and has given excellent results. Most of the earlier spectro- 
grams were taken on the Eastman 33 plate sensitized with pina- 
cyanol, or on the Ilford Special Rapid Panchromatic emulsion. 


* Mt. W. Contr., No. 432; Ap. J., 74, 188, 1931. 
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All the grating plates on which the D lines are visible were meas- 
ured for line displacements by Miss Burwell, and most of them were 
remeasured by her after the lapse of three months or more. On the 
remaining plates the second measure was by Mr. Merrill. The prob- 
able error of the mean displacement of the D lines derived from a 
single spectrogram appears to be about 2.2 km/sec for the 18-inch 
camera and 3.0 km/sec for the 1o-inch. 

Photometric tracings were made of all plates from which it ap- 
peared beforehand that a reliable result should be obtained. The 
principal material rejected may be classified as follows: (1) poor 
spectrograms; (2) lines of extremely low intensity; (3) early spectro- 
grams without photometric calibration; (4) dense exposures, poor 
for intensity measurements, although often excellent for displace- 
ments. As a result of this selection practically all the tracings were 
measured, and very few data were later rejected. Four-fifths of the 
spectrograms used for intensity measurements were taken with the 
18-inch camera, dispersion 34 A/mm. 

A set of standard circles for photometric calibration was impressed 
on each plate immediately after the star exposure, by means of a 
tube sensitometer having seven apertures of known area. The ex- 
posure in the sensitometer was in two parts: first for about 5 sec- 
onds; then, after the plate had been shifted half the diameter of one 
of the circles, a second exposure for about 12 seconds. This method? 
produced 21 standard densities in 3 overlapping groups of 7 points 
each, which define the calibration-curve with sufficient accuracy. 
The source of the calibrating light was a Mazda lamp shielded by an 
orange filter, which shuts out light of wave-length shorter than 
5535 A. Trials with a spectrum sensitometer showed that the III C 
emulsion yielded calibration-curves having sensibly the same slope 
from 5535 A to the effective limit of sensitivity of the plate at or be- 
yond 6400 A. 

The line intensities were determined in the usual way by measur- 
ing the widths of the lines at selected points below the level of the 
continuous spectrum and computing the equivalent widths by a 
simple process of numerical integration. Partially independent val- 
ues were obtained also from the measured central intensities by an 


2 Suggested by Mr. Wilson. 
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empirical method employed by M. Minnaert and B. van Assen- 
bergh for Fraunhofer lines, and by E. G. Williams‘ for interstellar 
K. For each line’ a smooth curve was drawn through the points ob- 
tained by plotting the observed central intensities against the equiv- 
alent widths found directly from the areas. The equivalent width 
corresponding to each measured central intensity was then read 
from the curve and combined with that obtained directly. 

A weight of from 1 to 3 was assigned each line at the time of meas- 
urement, in accordance with the observer’s judgment of the reliabil- 
ity of the resulting value. Since all the measuring and weighting 
were done by one person (Mr. Wilson), the system should be fairly 
homogeneous. 

One of the chief sources of error in measurements of equivalent 
line-widths usually lies in the difficulty of determining from the 
tracing the exact position of the continuous spectrum near the line. 
In the present investigation this difficulty was minimized by the use 
of the fine-grained III C emulsion and by the fact that few of the 
spectra are very narrow. Moreover, in these early-type spectra rela- 
tively little uncertainty is caused by neighboring lines. The presence 
of the alpha band due to oxygen in the earth’s atmosphere did, how- 
ever, cause some difficulty in the measurement of the interstellar line 
d 6284, especially on spectrograms taken at a considerable zenith 
distance. The violet wing of \ 6284 overlaps the head of the alpha 
band. Attempts to measure the total intensity of the blended fea- 
tures, and then to subtract a correction for the a head, determined 
as a function of zenith distance, did not yield satisfactory results. 
The method adopted was to measure the intensity of \ 6284 by it- 
self on the assumption that the blended violet wing has the same 
shape as the red wing. In a few stars with very weak interstellar 
lines the C 11 pair at 4 5889.97 and \ 5891.65 interferes with the meas- 
urement of D2; attention is called to several instances in the notes 
following the Catalogue (Table 6). 

In the measurements of the D lines we have, in general, assumed 
the contribution made by stellar D lines to be negligible. In types 


3 Zs. f. Phys., 53, 248, 1929. 
4 Mt. W. Contr., No. 487; Ap. J., 79, 280, 1934. 
5 D2 and D1 were taken together. 
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Ao and later, and in a few stars of earlier type, this assumption may 
not be precisely true; but we think the errors caused by the inclusion 
of stellar components are usually smaller than the accidental errors 
of observation, and of minor consequence. The question of the in- 
tensity of the stellar D lines in c stars, particularly those of classes 
B8 to A2, needs especial consideration. Most of these stars are so 
distant that the interstellar lines are strongly predominant and no 
direct measure can be obtained of the stellar lines. The indirect evi- 
dence from displacements and from correlations with other inter- 
stellar lines is that the stellar D lines are very weak. Some data are 
available for 8 Orionis cB8.and for a Cygni cA2, the brightest c stars 
of early type. The measured equivalent widths of D2 and D1 in B 
Orionis are 0.14 and 0.04 A, respectively. The D2 line is abnormally 
intense compared to D1, probably because blended with lines of 
Cu. If 8 Orionis is a typical ¢ star, its distance should be about 100 
parsecs; hence the interstellar D lines might be expected to have in- 
tensities of a few hundredths of an angstrom. In any event, the in- 
tensities of the stellar D lines are very low. In a Cygni the mean 
equivalent width of the D lines is about 0.5 A, but these lines appear 
to be chiefly interstellar. High-dispersion coudé spectrograms reveal 
clearly that the lines are composite, one component being deep and 
narrow and typically interstellar in appearance, and the other, shal- 
low and slightly diffuse, resembling the other stellar lines. The two 
components always overlap, but at times the relative displacement 
of the lines is sufficient to show that the stellar lines are not of great 
intensity. A rough estimate is that they do not contribute more than 
20 per cent of the total intensity. A more detailed investigation is 
desirable. 

From the results of nine plates of HD 224151 (Boss 6142), Table 1, 
it appears that the probable error of the measurement of a line of 
moderate intensity on a single good plate with a dispersion of 
34A/mm is about 8 per cent of the intensity for the D lines and about 
14 per cent for \ 5780 and \ 6284. For weak lines, as well as those on 
spectrograms below standard quality, errors of 20 per cent are prob- 
ably not uncommon. 

As a check on systematic error we have measured the intensities 
of the D lines on a few spectrograms of skylight. The appearance 
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of the solar D lines on our spectrograms does not differ very greatly 
from that of interstellar lines of comparable intensity. The standard 
method was used, of course, in the measurement of the control 
plates. The results from plates with the 18-inch camera, with which 


TABLE 1 
EQUIVALENT WIDTHS FROM INDIVIDUAL PLATES; HD 224151 
D2 Dr 5780 d 6284 
0.56A 0.44A 0.19 A 0.42A 
.64 .58 
.58 .70 
.60 . 26 .48 
.58 .46 .22 .50 
66 16 .48 
36 .26 .60 
62 .26 .54 
0.65 0.53 0.18 0.43 
Mean 0.62+.011} 0.50+.015 0.21+.011 0.52+ .020 
PE single plate + .033 + .044 O88 + .060 
TABLE 2 
MEASUREMENTS OF SOLAR D LINES 
Observer D2 Dr Remarks 
0.66A 0.46A High dispersion 
.58 44 High dispersion 
63 High dispersion 
Thackeray......... 65 A, High dispersion 
0.70 0.56 High dispersion 
ee 0.66 0.52 High dispersion 
Mt. Wilson........ _ 52 34 A/mm, 7 plates 
; Mt. Wilson........ 0.80 0.60 66 A/mm, 4 plates 


most of our stellar spectrograms were made, are in good accord with 
those obtained by observers using high dispersion, while those from 
the 10-inch camera plates are somewhat higher (Table 2). 

The only measurements of interstellar D lines with which we may 
compare ours are those of Beals® for the 15 stars listed in Table 3. 
The Mount Wilson values are systematically smaller by about 30 


5 M.N., 96, 661, 1936. 
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per cent, although the nearly linear grouping of the points in the cor- 
relation diagram (Fig. 1) indicates that in both series the accidental 
TABLE 3 
COMPARISON WITH MEASUREMENTS BY BEALS 
EQUIVALENT WIDTHS IN ANGSTROMS : 
BEALS Mount WILSON 
STAR 
D2 Dr D2 Dr 
0.34 0.29 0.24 0:22 
33656 0.82 0.61 0.50 0.38 
30861 0:35 0.26 0.42 0.24 
37043 0.19 0.14 0.14 0.07 
0.24 0.16 0.26 0.18 i 
1.08 0.86 O.77 0.57 
0.48 0.52 0.36 0.34 
A 
122 
4 
7 
0.8} 
fo) 4 e 
= ° 
4 
e 
° 
| 
e 
ree 
Cp 
; 
0.0 
0.0 0.4 0.8 1.2 1.6A 
Beals 


Fic. 1.—Equivalent widths of D lines measured at Mount Wilson and by Beals. | 
Dots D2, open circles D1. 


PLATE VI 


| 
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SPECTROGRAMS OF INTERSTELLAR LINES 
(a) HD 192163: mag. 7.4, sp. Ob; (b) HD 14143: 6.7, cB1; (c) x Aurigae, HD 36371: 
4.9, CB3; (d) w' Scorpii, HD 144470: 4.1, Bzn; (e) HD 21291: 4.4, cBo; (f) (g) 
HD 223385, Boss 6111: 5.6, cA2ea; (h) HD 193793: 6.8, Oa; (7) x? Orionis, HD 41117: 


4.7, cBrea; (7) HD 223987: 7.6, Bo. 
In (g) dots, left to right, indicate: \ 6278, head of the terrestrial alpha band; \ 6284; 


d 6347 Siu, stellar; \ 6371 Si 0, stellar. 
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errors are satisfactorily small. The source of the systematic error 
is not clear. The grating with which our spectrograms were taken is 
known to be particularly free from scattered light, while Beals’s spec- 
trograms were obtained by means of a prism. Hence it is not prob- 
able that scattered light can account for the difference in the re- 
sults. 


D3 De Di 


A5780 AS797 


Oz Atm. A6284 
4 


Fic. 2.—Tracings of interstellar lines in the yellow and red. (a) and () also show 
the stellar line D3. (a) w' Scorpii, HD 144470: mag. 4.1, sp. Ban; (6) HD 14143: 6.7, 
cBr; (c) (d) HD 223385, Boss 6111: 5.6, cA2ea. In (d) the feature marked O2 Atm is 
the head of the alpha band due to molecular oxygen in the earth’s atmosphere. Spectro- 
grams showing these lines are reproduced in Plate VI. 

OBSERVATIONS OF H AND K 

Most of the photographs of the H and K lines were made with a 
three-prism ultraviolet spectrograph having a camera of to inches 
focal length, dispersion at K, 24 A/mm. A 15-inch camera, disper- 
sion 17 A/mm, was employed for a few bright stars. A few early ob- 
servations were made with a one-prism spectrograph, dispersion 
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about 24 A/mm. More than three hundred spectrograms of two 
hundred stars were obtained, largely during the interval 1931-36. 
Most of the plates were taken with the aid of the 60-inch telescope; 
but for a number of the fainter stars, the 100-inch was employed. 
The Imperial Eclipse Soft emulsion HD 850 was used throughout 


K 


Fic. 3.—Tracings of interstellar H and K lines. (2) HD 193793: mag. 6.8, sp. Oa; 
(6) x? Orionis, HD 41117: 4.7, cBrea; (c) HD 223987: 7.6, Bo. Spectrograms showing 
these lines are reproduced in Plate VI. 


except for a few exposures on bright stars with the Eastman 33 or 
Eastman 40 emulsion. Most of the spectra were widened. 

The H line, as well as the K, was measured for radial velocity 
whenever possible, although the result was often given low weight 
because of the unfavorable influence of the neighboring He line. 
Nearly all the velocities depend on two measures of each plate by 


Mr. Sanford. 
The probable error of the velocity from a single plate, as deter- 


| H,He | 
b 

| 

| 
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mined from stars with numerous plates, appears to be somewhat less 
than 3 km/sec. For stars with two plates the mean arithmetic 
difference of the measured velocities is 4.5 km/sec. The mean arith- 
metic difference between a velocity determined at Victoria and the 
result from a single Mount Wilson plate of the same star is 4.2 km/ 
sec. The mean arithmetic difference between the first and second 
measure of the same plate for the first 50 entries in the Catalogue 
(Table 6) is 2.5 km/sec. These last three comparisons make it prob- 
able that 3 km/sec is a liberal estimate of the probable error of a 
single plate. 


TABLE 4 
WEIGHTS OF RADIAL VELOCITIES 
AVERAGE PE (Km/SEc) 
TypicaL NUMBER 
WEIGHT 
OF PLATES 
Ca Na 
1.8 3 or more good 
ere 3.0 2.2 1 good or 2 fair 
| 4.0 3.0 1 fair 


The intensity of the K line was measured by the same method as 
that employed for the D lines. Miss Grace Wilson made the tracings 
and drew the characteristic curves, while Mr. Sanford made the ac- 
tual measures. The standardizing exposures were impressed on each 
plate either with a tube sensitometer or by the spectroscopic method 
described by E. G. Williams.4 The accuracy of the determinations 
of equivalent width may be inferred from the fact that the mean 
arithmetic difference between two plates of the same star ranges 
from 11 per cent of the value for the strongest lines measured to 20 
per cent for the weakest. 


DESCRIPTION OF THE CATALOGUE 


The accuracy of the velocities recorded in the Catalogue (Table 6) 
varies so widely from star to star that it seems desirable to indicate 
the approximate weight of each entry. This has been done by use of 
the letters ‘“‘A”’ to “‘E,”’ whose significance is explained in Table 4. 


| 
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A number of values of the calcium velocity in column 8 of the Cat- 
alogue (Table 6), not accompanied by weight symbols, are measure- 
ments at other observatories, taken from Moore’s Catalogue of 
Radial Velocities. 

Similarly in the data for intensity, ‘‘A’’ indicates a weight greater 
than 4; “B,” weights 2 to 4; and “C,” weight 1, the weight assigned 
to a single spectrogram varying from 1 to 3. On those plates which, 
for various reasons,® were not run on the microphotometer, the line 
intensities were estimated by visual inspection under the micro- 


TABLE 5 


NUMBERS OF STARS OF VARIOUS SPECTRAL TYPES 
IN THE CATALOGUE, TABLE 6 


T No Additional 
Symbol 
53 17 18 
73 14 63 
22 21 16 
7 21 3 
215 76 109 


scope. These estimates do not lead to reliable numerical values but 
are useful for giving an approximate idea of the intensity, especially 
to indicate the spectra in which the lines are very weak. As nearly 
as can be specified, the equivalent widths for the K line or for o.5 
(D2 + D1),’ corresponding to the Roman numerals, are 
I..<o.04 A Line either not seen or a barely perceptible trace 

II....0.06+A Definitely present 

III....0.1 too.3 A Well marked 

IV..>03A Strong 


7 Pub. Lick Obs., 18, 1932. 
8 Mainly very low intensity of the lines. 


9 The estimates in the column headed ‘‘D2”’ are for the mean of the D lines. 


| 
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Except where otherwise indicated, stellar components are believed 
to add little to the observed intensities, and no corrections have been 
applied to the estimates. Thus an estimate of I or II means that both 
stellar and interstellar components are of low intensity. 

Of 400 stars in the Catalogue (Table 6), 254 are brighter than ap- 
parent magnitude 7.0, 124 are between 7.0 and 9.0, and 19 are fainter 
than 9.0. The distribution of types is indicated in Table 5s. 

Displacements of the sodium lines are recorded for 264 stars, meas- 
ured intensities for 213, and estimated intensities for 107. For the 
calcium lines the corresponding numbers are 174 (Mount Wilson 
values), 114, and 89, respectively. The intensity of \ 6284 is record- 
ed for 129 stars; of \ 5780, for 117 stars. 


| 
| 
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In Table 7 are collected all measurements of displacements of in- 
terstellar D lines made at other observatories that have come to our 
attention. The general agreement between these values and those 


determined at Mount Wilson is good. 


DISPLACEMENTS OF INTERSTELLAR D LINES MEASURED 


TABLE 7 


AT OTHER OBSERVATORIES 


HD Star Other Obs. Mt. W Ref. 
km/sec km/sec 
y Cas iI — 3B 2 
ere ¢ Per (+23)* +12A 2 
Ori +21C 2 
Ori +17 +20C 3 
e Ori +18 +20A 4 
Oni +17 +18B 4 
S Mon +11 +16A 3 
B Sco —9 —13B I 
A Cyg —17 —16B 2 


* Considered stellar by Miss Heger. 


t Miss Heger states that ‘the D lines unquestionably oscillate but there is not. suffi- 
cient evidence to show whether or not their amplitude is less than that of the other lines.” 


An bn 


Table 8 lists all stars for which the measured intensity of 0.5 
(D2 + Dr) or of \ 6284 exceeds 0.80 A, and in addition 7 bright 
stars selected because the detached lines in their spectra are un- 
usually intense for the apparent magnitude. These are strategic ob- 


. M. L. Heger, Lick Obs. Bull., 10, 59, 1919. 


. Ibid., p. 141, 1921. 


. J. S. Plaskett, Pub. Dom. Ap. Obs., 2, 287, 1924. 
. C.S. Beals, M.N., 96, 661, 1936. 
. J. A. Pearce, Pub. A.A.S., 8, 250, 1936. 


. N. T. Bobrovnikoff, Contr. Perkins Obs., No. 6, p. 29, 1937. 


jects to observe in a search for additional interstellar lines. 
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TABLE 8 
STARS WITH STRONG INTERSTELLAR LINES 
INTENSITY 
OTA MAG.| SPEC. 
1900 1900 
ge | 
2 
E.A. E.A E.A 
0"36™|+63°45'| 7.4/cB8 go"l+ 2°1...... 0.50 | 0.94 
37 | 8.8icAr 0.98 | 1.16 
44.6/+59 7 | QI |— 3 | 0.82 | 0.78 
ABAT 45.4/+63 14 | | 0.68 | 0.97 
OFO5%).... 1 24.6/+62 51 | 7.5\cBs(e)| 95 |+ 1 |..... | 0.89 | 1.00 
| 
12953.....| 2 1.7|\+57 57 | 5.9|cAzea |— 2]......| 0.85 | 0.65 
6.6/+58 6] |ror |— 2]......| 1.00] 0.81 
8.9|/-+57 50 | 7.8icAo |102 j— 2]...... 0.90 | 0.78 
8.9l+55 32 | 8.0] Bin |102 |— 4 ]|...... | 0.81 | 0.65 
9.8)/+56 34 | 7.2\cB2 102 |— 3] 0.72 0.96 | 0.60 
9.9/+56 36 | 6.4\cBrea |102 |— 3 | 0.75 0.82 | 0.62 
11.1/+63 58 7.0\cB8 0.87 | 0.85 
vik? 12.1|+56 40 | 6.7\cB3ea |102 |— 3 | 0.80] 0.87 | 0.58 
12.21+56 43 6. 7|cBr 102 |= 3 | 0:6) I.10 | 0.94 
14.8/-+56 47 | 6.5jcA2 |— 3]...... 0.81 | 0.88 
15.8/+56 47 | 7.5|cA2 |103 |— 3 ]...... ©.97 | 0.99 
15.9|+56 56 | 7.0\cB9 |103.|— 3 }.....- 0.98 | 0.78 
14947..... 19.5|+58 25 | 8.0] Osf |103 0.66 | 0.84 
22.8/+57 22 | 7.3\cA2 |r04 |— 2 ]...... 1.08 | 0.96 
15497... 24.61+57 15 | |104 |— 2 1.44 | 1.38 
27.21+-60 6 | 8.4] Bat |....-. 0.75 | 0.82 
| 36.3I+50 24 | 7.7\cCBo |+ 1 ]...... 0.85 | 0.88 
| | 3 8.21+56 45 | 5.9\cAo |109 |+ 1 |...... 0.61 | 0.87 
21389.....| 21.9/+58 32 | 4.8|cAoea |110 |+ 3 |..... 0.70 | 1.00 
5 50.9|+24 14 | 6.0\cBo |...... 0.46 | 0.92 
6 10.8/+23 46 | |+ 5 | 0.50] 0.54 | 0.88 
13.31+23 19 | 7.0] Bo |+ 5 |...... 0.47 | 0.86 
¢'Sco |16 47.0/—42 12 | 4.9/cBieq |311 0:73 | 
17 35.3|—33 27 | 6.7\cAgea |323 |— 3 |....-. 0.85 | 1.36 
2.6/—ar 28 | 2 |...... 0.63 | 0.84 
TOOASAS 19.6/—14 2] |345 |— 2] 0.60] 0.76 fore) 
183143.. 19 23.0/+18 5 | 6.9icBoea | 21 |— 1 |...... 0.68 88 
190603. . 20 0.7/+31 56 | | 37 | 0.39 0.68 | 1.02 
2.0/+35 37 | 9.5| Bne | 40 |+ 1] 0.40] 0.81 ]..... 
227836.....]...... 4.3\+35 50 |10.9| Bze | 41 [+ 1 |(0.46)| 1.43 |..... 
17.1|+43 32 | 6.8} Oa 49 |+ 3 | 0.49 | 0.70 | 0.84 
19.7|+40 26 | 7.0] Bo 46 |+ 1|0.78| 0.76 | 0.94 
207200. . vCep |21 42.6/+60 40 | 4.5\cA2 0.82 | 0.88 
45.6|+40 40 | 6.5\cAo 58 |—I0 ]...... 0.85 | 0.56 
22 19.3/+54 44 | 8.4/cB3 71 |— 2] 1.35 | 0.96 | 0.78 


* The intensities of the K line are from measurements by Beals, Williams, and Sanford. 
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TABLE 8—Continued 


INTENSITY 
HD STAR por Mac.| Spec. | 1 b 
me | 
2 

E.A E.A. E.A. 

23"44™0|+61°40'| 5.6;cA2ea | 83°} o°%}...... 1.27 | 1.15 
48.9/+60 18 | 7.0\cAoea 84 |— 1 ]...... 1.27 | 0.94 
49.2/+61 3] Bo 84 0.68 | 0.89 | 0.80 
49.7|+61 17 | 7.2!cB3ea | 84 ° .85 | 0.99 | 1.00 
| 58.8/+60 33 | 8.6} Bo 85 |— | 0.95 | 1.24 
58.9/+61 40 | 8.6) O8ea | 85 0.82 | 0.68 

Supplementary List of Bright Stars 

E.A. E.A. E.A. 

2905..... «Cas | 4.2)cBoea | 89°} 0°] 0.31 | 0.53 | 0.49 
3 21.0\+59 36 | 4.4/cCBg |+ 4 ]...... .48 
x Aur | 5 26.2/+32 7 | 4.9\cB3 1 .46 .40 
x?Ori | 5 58.0/+20 8 | 4.7\cBrea |157 |+ 1 .42 .4I .62 
@ |16 15.3|—25 21 3.11 Br (310° |...... 18 .41 
198478..... 55 Cyg |20 45.5/+45 45 | 4.9,cB2ea | 53 |+ 1 .50 42 .40 
23 58.3/+63 5 | 6.3\cB2ea | 85 |+ 1 | 0.67 | 0.64 | 0.64 
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THE SPECTRUM OF H, FROM } 4225 TO 2 4756 
NORTON A. KENT AND REGINALD G, LACOUNT 


ABSTRACT 

The wave-lengths of 435 lines of Hz have been measured on plates taken with a 
21-foot concave grating at the California Institute of Technology. The results have 
been compared with the measurements by Gale, Monk, and Lee, and by Finkelnburg. 

This paper is a continuation of data given in a previous paper.’ 
The apparatus employed was the same as that described therein, 
and the wave-lengths were determined in a similar manner. 

In Table 1, Nis the integer value of the wave-length in Angstroms; 
F, Finkelnburg’s fractional value; G, that obtained by Gale, Monk, 
and Lee; L, the Lacount? interferometer value; M,,, the writers’ 
weighted mean determination; V, the number of “‘sets’”’ of measure- 
ments entering into the calculation of M,,; PE, the probable error 
of the weighted mean, M,,; and v, the “vacuum”? wave-number’ in 
m~ obtained from M,, or from L, when the latter is present; } indi- 
cates that the line was broad; d, that the line was double and that 
the setting was made on the center of gravity; v, that there existed a 
weak, unresolved line on the violet side of the line measured. 

With the interferometrically determined wave-lengths, L, the 
grating determinations, M,, (shown in parentheses and based on the 
iron standards), are listed for comparison. The maximum deviation 
of M,, from L is 0.006 A, the average, 0.003 A. The Gale, Monk, and 
Lee interferometer determinations are marked by an asterisk in col- 
umn G of the Table. The average probable error of the mean for the 
values of M,, for the 435 lines listed is 0.001 A. 

The writers’ values agree, in general, more closely with those of 
Gale. The average difference between G and M,,, without regard to 
sign, is 0.008 A; that between F and M,, 0.014 A. 

A comparison of the 12 interferometer values, L, with M, reveals 

Ap. J., 84, 585, 1936. 

2 Obtained by Lacount with the same apparatus as that described in the article by 
Kent and Lacount, Phys. Rev., 51, 241, 1937. 

3 Calculated from Kayser’s table. 
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TABLE 1 


F 
26 790 
27 355 
28 860 
30 998 
32 683 
32 925 
33 408 
33. 893 
37 9°07 
42 476 
43 129 
43 348 
40 072 
46 997 
53 301 
54 040 
267 
62 794 
65 146 
67 go2 
77 475 
79 387 
80 724 
81 814 
86 068 
87 570 
89 637 
go 099 
- 425 
95. 667 
98. 034 
4303. 437 
849 
03 944 
05 577 
05 822 
06. . 273 
08. . 618 
II 715 
12 304 
12 924 
17 203 
18 O12 
18 261 
27 360 
27 939 


v 


2305891 


23605219 
2304864 
2304057 
2363304 
2361905 
2361784 
2301514 
23012608 
2350184 
2358923 
2350452 
23560085 
2355974 
23541006 
2353942 
2350450 
2350044 
23490394 
2349574 
2345533 
2345224 
2343910 
23423907 
2337173 
2330110 
2335381 
2334802 
2332501 
2331647 
2330543 
2330286 
2327403 
2327285 
2320007 
2323077 
2322829 
2322784 
2321982 
2321784 
2321538 
2320875 
2320269 
2318612 
2318244 
2317985 
2315605 
2315238 
2315117 
2310226 
2300924 


* Doublet just resolved in one case. Four others show an unresolved doublet with maximum den- 


sity at 0.921 (PE=0.002). 


: G I L Mw N PE . 
399 7 386 2 j 
670 682 2 002 
892 809 2 ool 
407 416 | 2 002 
818 824 3 002 
545 565 | 2 
125 139 4 002 i 
326 338 5 oor 
601 2 7O5 2 
988 o |......](47)002 
029 I 046 3 
778 790 2 003 
162 179 2 002 
480 474 4 003 
740 I ere 756 3 002 | 
816 I 818 2 002 
641 ge 643 2 
114 110 4 OOo! 
423 431 5 000 
423 429 6 000 
818 I ee 826 2 000 
276 3 282 (283 | 3 002) 
635 637 002 
701 716 5 
3904 O 401 2 
; 873 883 002 
004 002 | 4 002 4 
220 225 4 
362 368 2 002 
927] 5 | 934] (031 | 4 | 002) 
: 
’ 


THE SPECTRUM OF H, 313 
TABLE 1—Continued 
ny F G I L Mw N PE v 

628 619 620 6 oor 2307426 
532 519 524 | 6 2305880 
304 309 324 2 ool 2305455 
427 412 422 3 002 2296400 
ne 408 540 3 544 (sar | 2 oo1) | 2295808 
037 gI5 931 3 002 2291921 
002 990 (72)000 | 3 002 | 2286642 
307 403 402 4 2282778 
961 055 963 | 3 002 2282486 
194 181 IgI 3 002 2281325 
84.. 443 450 455 2280147 
86.. 201 249 260 3 000 2279208 
0905 084 og! 4 OO1 2277738 
go. 8090 goo 4 002 2276707 
738 726 742 4 oor 2276363 
101 093 100 | 3 002 2276178 
930 891 goo | 2 004 | 2275763 
ogI o81 093 3 ool 2273076 
242 2390 246 3 004 22720997 
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314 NORTON A. KENT AND REGINALD G. LACOUNT | 
TABLE 1—Continued 
r F G I L Mw N PE v 
487 478 496 3 oor 22606685 
650 618 628} 3 oor 2266617 
505 499 500 3 002 2265141 
209 218 217 5 002 2264773 
994 3 002 2264374 i 
485 480 477 2 003 2263615 
338 342 337 6 Oo! 2263174 
665 673 671 5 2262490 
490 404 403 5 000 2262070 
838 818 : 833 2 ool 2261896 
306 304 300 | 6 000 2261657 
928 gI5 920 2 OO! 2261340 
186 163 168 | 2 2261213 i 
goo 895 go2 5 oor 2258704 
738 750 749 | 3 002 2257852 
826 876 894 2 002 2254723 
278 278 271 2 oor 2254531 
217 206 205 2 2252533 
637 629 | 636] 3 002 2252314 
427 429 | 423 2 002 2251915 3 
875 852 857 | 3 004 | 2250174 
220 213 220 6 2249484 
245 246 246 5 oor 2248965 


Tt Measured with the interferometer. 


THE SPECTRUM OF H, 315 
TABLE 1—Continued 
aN F G I L Mw N PE v 

55° 353 004 | 2247797 
47. 928 932 O 931 4 2247607 
49. OO 307 2 000 2246882 
49. 512 502 514 | 2 004 | 2246808 
49. 919 6 oor 2246604 
50. 349 350 351 2 oor 2246385 
50. 817 816 828 | 3 2246144 
51. 357 3406 450} 3 002 2245882 
119 007 115 2 2245495 
503 482 407 2 2245303 
775 767 7761 3 oor 2245162 
169 141 158 6 002 2244970 
451 461 487 2 002 2243796 
664 666 686 3 003 2243696 
959 055 9057 2 002 2243559 
650 665 650 2 006 2243211 
872 851 rer 855 4 ool 2243107 
044 030 037 2 005 2243016 
568 577 577 2 002 2242744 
464 464 By 460 2 002 2241797 
189 190 177 4 002 2241437 
971 965f| 10 967 (963 | 4 000) | 2241042 
235 221 223 3 002 22390405 
209 280 285 2 oo! 2237369 
953 961 956 4 002 2235533 
269 261 258 | 5 002 2234383 
981 987 987 | 4 002 2232023 
075 083 | 2 003 | 2230481 


| 

5 

| 

? 

| | 

| 


316 NORTON A. KENT AND REGINALD G. LACOUNT 
TABLE 1—Continued 
r F G I L Mu N PE v 

873 832 842 2 2228613 
145 103 109 Be 2227983 
824 813 816 8 oor 2227632 
574 548 549 2 000 2227269 
214 234 248 2 002 2226922 
458 9 453 (450 2 000) | 2226325 
741 738 736 2 003 2225688 
693 688 692 7 002 2224720 
310 313 323 2 002 2223013 
666 672 679 3 003 2223242 
096 IOI 114 2 000 2223027 
569 577 580 3 003 2222707 
108 10 |. 56... 106 7 ool 2222537 
548 554 555 4 002 2221821 

046 050 053 2 003 2221575 
595 597 602 2 000 2220811 
g6o g60 959 4 oo! 2220635 
598 600 614 2 000 2220311 
4604 407 | 476 3 002 2219394 
928 g16 O27" | 002 2219172 
635 631 631 4 003 2218825 
076 103 III 2 000 2217112 
898 904 905 3 005 2216230 
705 690 Be 696 4 004 2215842 
837 828 828 2 002 2214795 
166 182 192 2 000 2214120 
579 562 568 | 4 002 2213942 
887 898 904 2 ool 22127097 
969 959 965 2 ool 2211788 
323 314 | 316 2 002 2210638 
184 193 195 3 oor 2210209 
149 139 7 145 (139 ) | 2209745 


THE SPECTRUM OF #H, 317 
TABLE 1—Continued 
r F G I L Mw N PE v 

085 079 Cha 076 3 003 2207339 
240 266 272 2 000 2207243 
clade 193 198 2 002 2206305 
053 047 Be 950 2 ool 2205936 
151 157 158 4 ool 2204865 
620 627 6 627] 4 2204637 
go2 goo 2 904 2 2204016 
880 899 2 003 2203048 
114 118 121 2 oo! 2201484 
928 954 I 959 2200594 
702 692 691 4 ool 2200239 
940 967 966 2 000 2198171 
goo 896 806 5 oor 2197237 
424 423 426 3 003 2196017 
164 |. 160 | 3 000 | 2195181 
125 125 135 | 4 002 2193748 
192 206 217 2 oo! 2192266 
732 728 736} 4 Oo! 2190575 
759 753 778 2 002 2188637 
setae 133 131 | 6 oor | 2188467 
616 617 620 | 3 2188233 
235 225 229 | 5 2186506 
880 878 880 | 4 2184762 
023 O15 o16 5 oor 2183742 


318 NORTON A. KENT AND REGINALD G. LACOUNT 


TABLE 1—Continued 


r F G I L Mw N PE v 
462 454 458 | 4 002 2183055 
, ee 996 904t| 10 | 904 (9903 | 4 | 003) | 2182799 
546 538 539 | 4 oor 2182063 
596 592T| 10 504 (504 | 4 002) | 2181561 
680 678 680 | 4 002 2178668 
581 579 596 2 003 2176336 
216 206 209 4 ool 2174626 
495 493 496 | 4 2174017 
417 403 409 | 4 2169811 
580 588 591 4 003 2166434 } 
524 528 528 | 4 oor 2165056 
304 304 | 310 | 5 2164690 
6590 600 598 | 3 2164554 
315 311 tf3190 | 5 002 | 2161410 
584 \s03 | 5 002 | 2161381 
984 |... - 988 | 4 002 | 2160163 
732 708 716 2 002 2159823 
703 673 BH 689 2 000 2159369 
069 074 103 2 003 2158710 
847 849 850 |} 5 002 2158362 
oS ere 030 032T| 10 031 (035 | 3 002) | 2157346 
891 898 895 | 2 2154020 


tf Faint broad line bridges region between. 
§ Broad, and interferes somewhat with 4634.035. 


THE SPECTRUM OF H, 319 
TABLE 1—Continued 
ny F G I L Mw N PE v 
007 (53)007 | 5 oor | 2148548 
188 196 7 204 | 4 2148457 
064 056 057 3 002 2148064 
861 861 861 2 000 2146309 
790 784 790 2 002 2145881 
te 390 305 402 5 002 2145139 
391 402 402 4 002 2144679 
812 811f| 8 810 (816 | 2 003) | 2144031 
595 585 Oo 587 4 2142755 
796 791 801 4 002 2141739 
976 | (70)o15 (7o)o20 | 2 004 | 2140721 
7). 847 808 811 2 000 2140359 
301 30st, 7 308 (309 | 2 003) | 2140131 
710 710 718 2 000 2139943 
093 0g2 093 4 002 2136570 
354 354 364 | 3 2135534 
341 341 348 | 5 2135085 
650 654 656 4 oor 2134033 
145 142 143 4 002 2133355 
774 772 784 | 4 003 2133005 
439 438 442 4 002 2132310 
040 040 048 oor 2130671 
723 723 729 3 000 2126284 
581 562 574 4 002 2125902 
502 489 500 | 2 002 2123226 
535 536f| 10 538 (538 | 2 004) | 2122759 
064 067 064 3 oor 2122071 
Q22 930 9290 | 4 002 2120781 
357 308 330 2 ool 2120601 
836 810 838 | 3 002 2118575 
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TABLE 1—Continued 


d F G I L Mw N PE v 
039 043T| 10 |...... 043 | 5 oor | 2118483 
003 000 003 2 000 2117603 
539 542 548 | 3 003 2117359 
O4I 30. 035 | 4 002 2116692 
BA 850 838 847 | 3 000 2116328 
828 820 827 4 002 2115889 
759 754 764 3 002 2113234 
992 9857 985 | 4 002 2108679 
094 109 114 4 002 2108176 
al 780 780 781 5 002 2107880 
391 383 385 4 002 2107612 
323 314 273 2 005 21060755 
042 968 (45)077 | 2 007 2106460 
581 575 578 2 003 2103978 
376 405 403 2 4 002 2102285 


that the sum of the positive values of L — M,, is 0.020 A and that 
the sum of the negative is 0.012 A, showing an average deviation of 
0.003 and 0.002 A, respectively, which indicates an approximately 
random distribution of differences. For L—G the sums are +0.034 
and —o.0o02 A; for L — F, +0.066 and —0.024 A. 

It is the writers’ intention to publish the wave-lengths of lines in 
other regions as soon as possible. 


Grateful acknowledgment is made to Dr. G. H. Dieke, who, while 
in California, assisted the senior writer materially in the long ex- 
posures necessary; to Dr. W. Uyterhoeven for a brief continuation of 
the work after the senior writer’s departure; to numerous Boston 
University students who have assisted in the computations; to the 
Rumford Committee of the American Academy of Arts and Sciences 
for its generous financial assistance; and to the Elizabeth Thompson 
Science Fund for the loan of a calculating machine. 


CALIFORNIA INSTITUTE OF TECHNOLOGY, 1927 
AND 
Boston UNIVERSITY 
June 1937 
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PHOTOCHEMISTRY OF PLANETARY ATMOSPHERES 
RUPERT WILDT 


ABSTRACT 

All polyatomic molecules detected in planetary atmospheres are highly sensitive to 
ultraviolet solar radiation. Their photochemical decomposition must be followed by 
secondary chemical reactions, reuniting the products of dissociation, in order to main- 
tain the observed stationary composition of the atmospheres. The stability of methane 
and carbon dioxide is reasonably accounted for, but that of ammonia remains unex- 
plained at present. The existence of methane and the failure to detect any other hydro- 
carbon indicate the presence of plenty of hydrogen on the giant planets, and most 
probably there is no carbon monoxide, carbon dioxide, or hydrogen sulphide at their 
surfaces. Small amounts of sundry molecules, not observed so far, are certain to be 
formed in consequence of the photochemical processes, and eventually these compounds 
may be detected under favorable circumstances. 

In a short note attention was called to the role photochemical 
processes play in planetary atmospheres, but information about the 
very reactions involving the main atmospheric constituents was still 
scanty then.’ A more detailed discussion seems profitable now, since 
numerous investigations have been completed meanwhile, dealing 
not only with the absorption spectra of polyatomic molecules in the 
far ultraviolet and the corresponding photochemical dissociations 
but also with the secondary chemical reactions of free atoms and 
simple radicals. These secondary reactions are of great importance 
to an understanding of the net result the primary photochemical 
process leads to. They are difficult to predict, and the final elabora- 
tion of a net reaction scheme, fitting the experimental results, is an 
intricate problem of chemical kinetics, still controversial in many 
cases. 

The principal task, from the astrophysical point of view, is to 
explain how a stationary state can be maintained in atmospheres 
containing molecules highly susceptible to photochemical decomposi- 
tion by ultraviolet solar radiation, like ozone, carbon dioxide, am- 
monia, methane, etc. For the ozone in the earth’s stratosphere this 
problem may be regarded as solved quantitatively, at least with re- 
spect to the order of magnitude of the stationary concentration of 


ozone. For the other gases mentioned and for some related com- 


™R. Wildt, Nature, 134, 418, 1934. 
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pounds, apparently less prominent, a qualitative solution will be 
attempted here. Which primary photochemical process is prepon- 
derant in a mixture of atmospheric gases depends on the energy 
distribution of the incident radiation and on the respective absorp- 
tion coefficients, and, therefore, indirectly on the atmospheric level. 
For example, the decomposition of ozone by the absorption of visual 
solar radiation (Chappuis-bands) shifts the photochemically main- 
tained stationary state of the stratosphere to lower ozone concen- 
trations, compared to the hypothetical state produced by the action 
of ultraviolet sunlight alone. On the other hand, carbon dioxide, oc- 
curring in traces in the troposphere, is prevented from decomposing 
photochemically by the protective ultraviolet absorption of the vast 
mass of oxygen in the upper atmosphere. Evidently, careful con- 
sideration is to be given to the relative abundance of the atmospheric 
constituents and to their vertical distribution. The variety of sec- 
ondary reactions increases greatly with the number of possible part- 
ners. For this reason the proposed reaction schemes leave room for 
some criticism. But the presence of atmospheric constituents not 
detectable by spectroscopic means may be postulated or refuted by 
the ensuing secondary reactions, i.e., by purely chemical arguments. 
It is indeed believed that the presence of a considerable amount of 
hydrogen in the atmospheres of the giant planets can be established 
in this way, which fortunately corroborates certain general cosmo- 
gonic conceptions. 

In order to illustrate the photochemical efficiency of a continuous 
spectrum like that of sunlight, it is appropriate to express the energy 
density as a function of the frequency of radiation. The familiar 
form of the emission power of a black body 


has to be replaced by 
h hy —I 
E, = 2x (ar (2) 


It is easy to derive the energy maximum on the frequency scale 


(3) 


| 
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where a = 2.8214 is the root of the transcendental equation 
a 
= (4) 
a 


On introducing the customary wave-numbers (unit cm™), the nu- 
merical value of the constant is a* = 1.97 cm™ degree, which re- 
veals that the maximum number of quanta of solar radiation (T = 
6000° K) is close to 11,820 cm™, or 8450 A. By substituting 


= xand—’- =y, (5) 


This normalized function is displayed by Figure 1, the abscissae be- 
ing proportional to the wave-numbers and the ordinates to the num- 
ber of quanta. Whether the energy distribution of solar radiation 
follows a black-body curve throughout the far ultraviolet cannot 
be ascertained, and the existence of an ultraviolet ‘“‘window,”’ 
situated beyond the Lyman continuum, has been much discussed 
recently. Nevertheless, Figure 1 illustrates two facts. On a fre- 
quency scale the Planck function has its steeper gradient on the in- 
frared side of the energy maximum, and the number of quanta 
decreases rapidly throughout the range from 40,000 to 80,000 cm™ 
for solar radiation. 

To start with the gaseous elements, there is no reason to suppose 
that their photochemical dissociation makes any significant con- 
tribution to the chemistry of planetary atmospheres; oxygen is to 
be excepted from this statement and will be dealt with later. The 
existence of free halogen molecules can be ruled out for their great 
chemical affinity. The monatomic noble gases have very high ion- 
ization potentials, and even in the ionized state they seem to be 
chemically inert. The limiting wave-lengths of the photochemical 
decomposition of molecular nitrogen and hydrogen are 1370 A and 
850 A, respectively. Since all planetary atmospheres are known to 


hy =2 
kT 
Planck’s function (2) may be normalized: 
} 
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contain polyatomic molecules strongly absorbing much greater 
wave-lengths, both hydrogen and nitrogen are probably well shield- 
ed. The best chance for the photochemical decomposition, or at 
least excitation, of molecular nitrogen seems to be given in the at- 
mosphere of Mars, which should contain the smallest amount of 
polyatomic molecules among all planetary atmospheres, judged by 


[ 


oW 


so 


40 ae 40 «0 so 70 @o go 400 


Fic. 1.—Graph of the normalized Planck function, eq. (6). The upper and the 
right-hand scales refer to the almost straight curve, giving log'® y as function of x in the 
range from x = 10.0to x = 20.0. Fora temperature of 6000° K, the following abscissae 
and wave-numbers correspond to one another: 


x = 2.39 N = 10,000 cm™! X= 9.45 N = 40,000cm™! 
77 20,000 14.32 60,000 
7.16 30,000 Ig.00 80,000 


the spectroscopic evidence. Of course, all planets are likely to be 
surrounded by ionospheres similar to the terrestrial one, but their 
masses should be exceedingly small compared with the bulk of the 
atmospheres. The consideration of ionospheric phenomena on other 
planets would go beyond the scope of this paper. But it seems worth 
while to point out that luminescence of an aurora-like character may 
also occur in the atmospheres of other planets. In the experiments 
made by Terenin and his collaborators,’ the photochemical decom- 
position of various inorganic and organic compounds produces ex- 


2H. Neuimin and A. Terenin, Acta physicochemica U.S.S.R., 5, 465, 1936. 
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cited radicals, like OH, CN, NH, which exhibit strong fluorescence 
both in the ultraviolet and the visual region.’ It is very tempting 
to ascribe to atmospheric luminescence the discrepancies between 
the planetary diameters observed in the light of different wave- 
lengths. But, as long as the physiological and photographic influ- 
ences affecting the planetary diameters have not been cleared up 
completely, every suggestion of this kind seems premature. 

The photochemical decomposition of oxygen molecules yields un- 
excited atoms alone or 50 per cent of excited atoms, according to 
the wave-length employed: 


< 1950 A: O, + hv O + O (3P), 
< 1760 A: O, + hv O BP) + O0('D). 


In the earth’s stratosphere this process is followed by partial recom- 
bination and by the formation of ozone, which subsequently decom- 
poses, both photochemically and thermally. The terrestrial ozone 
equilibrium may be dismissed by referring to the most recent treat- 
ment along theoretical lines.4 No other planetary atmosphere is 
known, thus far, to contain spectroscopically observable amounts of 
either oxygen or ozone. But it is of interest to note that in an at- 
mosphere containing only very little oxygen, say about 5 g cm”, the 
ozone layer would be formed close to the ground and, therefore, 
might attack chemically the planetary crust. Such may have been 
the case on Mars, in an earlier stage of evolution. For this and other 
cosmochemical implications, see a paper giving a more detailed dis- 
cussion of oxygen and ozone in planetary atmospheres. 


PHOTOCHEMISTRY OF SIMPLE HYDROCARBONS 


The low-pressure bands of methane start at \ 1300, and at higher 
pressure they extend to a limit around \ 1450. From there to \ 850 
the absorption is entirely continuous, and therefore all excited states 
up to the first ionization potential (approximately 15 volts) appar- 

3 Such phenomena presumably have some bearing upon the luminescence of comet- 
ary gases. The writer intends to take up the chemistry of cometary matter in a later 


paper. 
40.R. Wulf and L. S. Deming, Terr. Magn. and Atm. El., 41, 299 and 375, 1936. 
SR. Wildt, Nachrichten Gesell. Wiss. Gottingen (Math. Phys. KI. [N.F.]), 1, 1, 1934. 
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ently are repulsive. The low-pressure absorption of ethane con- 
sists of extremely diffuse bands starting at \ 1350 and going down 
below \ 1000. With increase in pressure the limit moves to \ 1500 
or even to \ 1600. On account of the extreme difiuseness of the 
ethane bands it is assumed that, as in methane, all upper electronic 
states are unstable.’ In ethylene, absorption has been observed 
from about A 1980 to A 1200, which consists of numerous well-defined 
bands.” * For acetylene, rather weak absorption has been reported 
between A 3150 and X 2870, and strong photochemical action is 
known to result from absorption in the region between \ 2300 and 
1900, where the bands are perfectly sharp.7° A third extensive 
system of more or less diffuse bands covers the range from A 1520 to 
A 1050. Naturally, the observed limits of absorption depend very 
much on the gas pressure and the length of the absorbing layer. Ac- 
cording to another investigator,’ the continuous absorption of 15 cm 
of methane under standard conditions extends up to A 1820. In any 
case, the sequence CH, — C.H¢ — C.H, — C.H, indicates the shift 
of the absorption limit to greater wave-lengths. Beside methane no 
other hydrocarbon has thus far been identified by the spectroscope 
in the atmospheres of the giant planets. But ammonia has been 
found on both Jupiter and Saturn. Since this gas absorbs from A 2300 
downward and, therefore, partly protects the methane, the following 
discussion of the decomposition of methane has direct bearing only 
upon the atmospheres of Uranus and Neptune, where all ammonia 
is probably frozen out because of the lower temperature. The com- 
plications arising from the fact that the products of the simultaneous 
decomposition of methane and ammonia are likely to react on each 
other will be taken up in a later section. Adel and Slipher™ estimate 
the amounts of methane on Uranus and Neptune to be approximate- 
ly 4and 25 mile-atmospheres, which may be a lower limit only for the 
possible pressure-broadening of the comparison spectra. By division 


6 A. B. F. Duncan and J. P. Howe, J. Chem. Phys., 2, 851, 1934. 
7W. C. Price, Phys. Rev. (2), 47, 444, 1935. 

8 G. Scheibe and H. Grieneisen, Zs. f. phys. Chem., B, 25, 52, 1934. 
9G. B. Kistiakowski, Phys. Rev. (2), 37, 277, 1931. 

10 A. Rose, Zs. f. Phys., 81, 758, 1933. 

™ Phys. Rev. (2), 47, 787, 1935. 
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of these numbers through the height of the homogeneous atmosphere 
(RT/gM) of methane one gets the partial pressure of methane at the 
base of the atmospheres, provided these consist exclusively of 
methane. If, however, as is most likely, the atmospheres of Uranus 
and Neptune contain a large surplus of hydrogen, this gas would de- 
termine the effective distribution of atmospheric density, and the 
height of the homogeneous atmosphere of hydrogen would be the 
appropriate divisor (about 30 miles for Uranus and 25 miles for 
Neptune). On this assumption the partial pressure of methane may 
be o.1 and 1 atmosphere for Uranus and Neptune, respectively. For 
Uranus this partial pressure is of the order of the vapor pressure of 
methane at the generally assumed surface temperature of the planet 
(go° K). For Neptune the surface temperature ought to be at least 
110° K in order to maintain the observed amount of methane in the 
gaseous state. This surface temperature seems by no means impossi- 
ble, in view of the powerful infrared absorption of 25 mile-atmos- 
pheres of methane. 

On irradiation of methane by Schumann ultraviolet light, the 
primary process is the formation of hydrogen atoms and methyl 
radicals; 


CH,+ hv H+ CH,. 


The detailed investigation of the quantum yield and of the final re- 
action products is still under way, but a preliminary announcement 
has been made which states that in pure methane the secondary re- 
actions produce mainly H, and C.H,, also small amounts of C.H¢6 
and traces of C,H, and of higher gaseous hydrocarbons.” Most for- 
tunately, sufficient information is available about the behavior of 
both hydrogen atoms and methyl radicals which have been produced 
and studied independently from the photochemical primary process. 
These experiments lead to the expectation that the foremost second- 
ary reactions should be recombinations by way of triple collisions: 


H+H4+M-—H,4+M, 
CH,+ CH,+M 


"2 W. Growth and H. Landenklos, Die Naturwissenschaften, 16, 796 and 828, 1936. 
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provided the reactions go on in an inert atmosphere, which does not 
furnish any chemical acceptors to the radicals. Neither hydrogen 
atoms nor free methyl groups react with surplus methane, according 
to Bonhoeffer and Harteck.'’ The recombination of free methyl 
groups, produced in an atmosphere of helium or nitrogen, has been 
followed quantitatively by Paneth.'* When he replaced the inert 
carrier gases by hydrogen, the formation of ethane was suppressed 
almost completely, and methane was re-formed by the reaction 


ON, +H, +2. 


For further data concerning this reaction, consult a recent paper by 
Patat.'® It would seem now that for maintaining a stationary con- 
centration of methane the presence of an ample supply of hydrogen 
would be a sufficient guaranty. This has indeed been concluded both 
from general cosmogonic reasoning” and from more specific evidence 
based upon the interpretation of the density distribution in the in- 
terior of the giant planets.'? Therefore it may be held that the pure- 
ly thermal reaction between CH, and H, frustrates the photochemi- 
cal decomposition of methane in the atmospheres of Uranus and 
Neptune. 

This statement is not invalidated by the observations, quoted 
above, on the formation of higher hydrocarbons during the photo- 
chemical decomposition of pure methane. This may result from sev- 
eral processes of the following type, sometimes called “‘dispropor- 
tionation”’: 


and from subsequent additions, like 


C+, «CM, 


3K. F. Bonhoeffer and P. Harteck, Grundlagen der Photochemie, 1933; K. H. Geib 
and P. Harteck, Zs. f. Phys. Chem., (Ser. A), 170, 1, 1934. 

™ F, A. Paneth, W. Hofeditz, and A. Wunsch, J. Chem. Soc., 1, 372, 1935. 

SF. Patat, Zs. f. phys. Chem., B, 32, 274, 1936. 

%*D. H. Menzel, Pub. A. S. P., 42, 228, 1930. 

7R. Wildt, Nachrichten Gesell. Wiss. Gottingen (Math. Phys. KI. [N.F.]), 1, 67, 
1934. 
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Also it should be noted that these higher hydrocarbons are likely to 
decompose photochemically in their turn, all the more since they ab- 
sorb at greater wave-lengths than methane. Acetylene polymerizes 
under the influence of ultraviolet light, and the reaction product is 
precipitated in solid form. Therefore, the absence of acetylene from 
planetary atmospheres is easily accounted for. To the writer’s 
knowledge, the photochemical decomposition of ethane has not yet 
been studied, but it may be presumed that a breaking-up of the C-C 
bond would be found to be the primary process. Then the methyl 
radicals would react with surplus hydrogen, forming methane, and 
ethane would prove to be photochemically unstable in a hydrogen 
atmosphere. Much work has been done on the photochemical de- 
composition of ethylene," but the interpretation of the experimental 
results is still beset with difficulties. Pending evidence to the con- 
trary, it would not seem likely in any way that ethylene could resist 
the combined action of hydrogen and ultraviolet radiation. More- 
over, another type of secondary reaction will tend toward the de- 
struction of all higher hydrocarbons in the planetary atmospheres. 
As pointed out before, hydrogen atoms originate in the re-formation 
of methane and disappear by recombination, which is a comparative- 
ly slow process for the triple collisions involved. Therefore, a small 
stationary concentration of hydrogen atoms ought to exist in the at- 
mospheres of the giant planets. But hydrogen atoms are known to 
attack violently all hydrocarbons, methane excepted, and to cause 
hydrogenation, breaking of carbon chains, and, principally, dehy- 
drogenation.'® Even if this process be a very slow one, it should have 
led to complete destruction of the higher hydrocarbons within a 
time which is short compared with the age of the planets. So it is 
very gratifying to realize that methane holds a rather exceptional 
position with respect to its stability. Otherwise one would fail to 
understand the absence from planetary atmospheres of other simple 
hydrocarbons, all of which have vapor tensions only slightly smaller 
than that of methane and, therefore, could not possibly be frozen out 
completely. The formation of methane in the cooling of a primeval 

*®R. D. McDonald and R. G. W. Norrish, Proc. R. Soc., London, A, 157, 480, 1936. 


9 K. F. Bonhoeffer and P. Harteck, Zs. f. phys. Chem., 139, 64, 1928; H. von Warten- 
berg and G. Schultze, ibid., B, 2, 1, 1929. 
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atmosphere rich in hydrogen and carbon oxides, as suggested by 
Russell,?? implies the simultaneous synthesis of small amounts of 
higher hydrocarbons, and apparently it is only by the photochemical 
processes and secondary reactions outlined above that the other 
hydrocarbons have been destroyed during a more advanced stage of 
evolution. It may be added here, too, that hydrogen sulphide, in 
addition to its photochemical decomposition, is rapidly disinte- 
grated by hydrogen atoms, solid sulphur being precipitated. 


PHOTOCHEMISTRY OF AMMONIA AND RELATED COMPOUNDS 


It is well known that the absorption spectrum of ammonia at ordi- 
nary temperatures and pressures starts at about \ 2200 and extends 
into the far ultraviolet. With pressures of the order of 1 mm or less, 
the bands down to X 1665 are all diffuse because of predissociation. 
Below A 1665 all the bands are very sharp and show rotational fine 
structure. True continuous absorption, in distinction to that pro- 
duced by pressure broadening, does not begin until about Xd 1200. 
At \ 1150 and below, the continuous absorption is so strong that no 
more bands could be measured accurately. Sharp bands exist, how- 
ever, as far down as A 1085.7" Two ways of photochemical decom- 
position have been observed in the laboratory, and a third one seems 
to be energetically possible, namely: 


\ < 2200A: NH; + hv H + NH., 
4 <1700A: NH, + NH,*, 
< 1900 A: NH, + hv H,+ NH*. 


The NH,* are excited radicals which show strong fluorescence be- 
tween \ 4000 and X 6000, emitting the so-called a-bands of ammonia 
observed in discharges through streaming ammonia vapor. The 
resonance bands \ 3360 and A 3370 of the molecule VH have been 
sought in vain. So the third process, which would cause the emission 
of the NH bands, can be ruled out because of experimental evidence. 
The first process should be preponderant in the atmospheres of 
Jupiter and Saturn, where the second is probably inhibited by the 


20 Science, 81, 1, 1935. 
2 A.B. F. Duncan, Phys. Rev. (2), 47, 822, 1935. 
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competing absorption of methane in the very range of wave-lengths 
required. Now the quantum yield of the first process is about 0.25 
at room temperature and varies only slightly with the partial pres- 
sure of ammonia. Evidently ammonia is re-formed at a considerable 
rate, by the secondary reactions. But it has not yet been established 
beyond doubt which course the secondary reactions take, despite ex- 
tensive experimental investigations. In all experiments part of the 
ammonia is destroyed irreversibly, the final products being molec- 
ular hydrogen and nitrogen. No reaction scheme can be conceived 
at present which would prevent any irreversible decomposition of 
ammonia and, therefore, would serve to maintain the presumably 
stationary composition of the atmospheres of Jupiter and Saturn. 
Even the presence of surplus hydrogen is of no direct help. By anal- 
ogy to the re-formation of methane, one might expect a reaction 


NH, + H, > NH; +4, 


which would imply a dependence of the quantum yield on the par- 
tial pressure of hydrogen admixed to the ammonia. The early ex- 
periments by Warburg prove the contrary. With due regard to the 
entire experimental material, an opinion may be ventured here as to 
the direction in which a solution of the astrophysical problem may be 
sought. But it is realized that any final verdict needs the support of 
further experiments. 

Farkas and Harteck” succeeded in measuring the stationary con- 
centration of hydrogen atoms in decomposing pure ammonia and 
explained their results by the following set of secondary reactions: 


a) NH,+H=NH,, 
NH, 
b NH, + NH; 
Ha+NH, 
c) 


At room temperature the equilibrium (a) is entirely on the side of 
the ammonium radical. The quantum yield is determined by the 
ratio of the alternative reactions (b) and (c), and in order to get the 


22 Zs. f- phys. Chem., B, 25, 257, 1934- 
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observed value (0.25) it has to be assumed that (6) proceeds three 
times faster than (c) at room temperature. Melville? made the in- 
teresting observation that in a mixture of ammonia and hydrogen 
the quantum yield of the photolysis of ammonia is very efficiently 
reduced by the addition of hydrogen atoms. Farkas and Harteck 
suggest that this effect is caused by the fact that in Melville’s ex- 
periments NH, and NH, are no longer present in approximately 
equimolal amounts, on account of the increased association of VH, 
with the added hydrogen atoms, and that, therefore, the surplus 
NH, may enter reactions like the following, or similar ones: 


NH, + NH NH, + N&,, 


consuming the fatal NH molecules which otherwise would decom- 
pose irreversibly. Now it will be remembered that a stationary con- 
centration of hydrogen atoms is set up in a mixture of hydrogen and 
methane on irradiation by ultraviolet light. Thus it may be that the 
very presence of methane in the atmospheres of Jupiter and Saturn 
prevents the ammonia from decomposing irreversibly, by a mecha- 
nism closely related to that which operates in Melville’s effect. 
Whether this is true or not is hard to ascertain at present for lack of 
knowledge of the extinction coefficients of methane and ammonia in 
the overlapping region of absorption, among other reasons. 

It should not be overlooked that other investigators** have pro- 
posed a rather different reaction scheme, involving the intermediate 
formation of hydrazine molecules. Since they did not measure the 
stationary concentration of hydrogen atoms during decomposition, 
their interpretation is necessarily at variance with that given by 
Bonhoeffer and Harteck and can scarcely claim the same significance. 
Hydrazine is decomposed by ultraviolet light too: 


< 2500A: N.H,+ hv NH,+ 
\ < 1600 A: N.H,+ hv > NH, + NH,*, 


where the excited radical emits the a-bands of ammonia referred to 
before. The so-called Schuster bands of ammonia, \ 5643 and A 5681, 
23 Trans. Faraday Soc., 28, 885, 1932. 
24R. A. Ogg, P. A. Leighton, and F. W. Bergstrom, J. Amer. Chem. Soc., 56, 322, 
i934. 
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are now ascribed to a transition from an excited stable state of the 
hydrazine molecule to a repulsive one.?5 Their absence from the 
fluorescence accompanying the photochemical decomposition of hy- 
drazine by Schumann ultraviolet light, as established by Terenin, 
does not favor the conception of a transitory formation of hydrazine. 
Besides, only rapid destruction by a subsequent reaction could pro- 
tect hydrazine from being frozen out. The exceedingly low-vapor 
tension of this compound at the surface temperatures of Jupiter and 
Saturn is by far the strongest argument against its appearance in the 
foregoing scheme of the photolysis of ammonia. However, the oc- 
currence of methylamine in the atmospheres of the same planets 
should be envisaged. Traces of the radicals CH, and NH, will cer- 
tainly recombine, and since in the range of temperature with which 
we are concerned here the vapor tension of CH,- NH, is about one- 
tenth that of ammonia, there is little danger of condensation. 
Whether these traces can be detected by the spectroscope is a differ- 
ent question. It would be worth while to look for the rotation-oscil- 
lation bands of methylamine in the spectra of Jupiter and Saturn, 
although they may conceivably be difficult to disentangle from the 
heavy absorption of methane and ammonia. The electronic absorp- 
tion spectrum of methylamine extends from \ 2400 downward, but 
the corresponding photochemical process has not yet been investi- 
gated. 

Furthermore, researches should be noted here concerning the re-. 
action of decomposing ammonia with ethylene or carbon monoxide. 
Although these experiments are of an orienting character only, most 
probably they preclude the existence of either of these gases in the 
atmospheres of Jupiter and Saturn. The hydrogen atoms originating 
in the photolysis of ammonia are found to induce polymerization of 
ethylene.” In principle this is the same process as discussed before, in 
view of the hydrogen atoms formed by the reaction of CH, with H.,. 
The photochemical reaction of ammonia with carbon monoxide’? 
yields formamide (HCO-NH,) as a primary product, which then 
undergoes secondary reactions and decomposes photochemically too. 


2*R. W. Lunt, J. E. Mills, and FE. C. W. Smith, Trans. Faraday Soc., 31, 792, 1935. 
*H.S. Taylor and H. J. Emeléus, J. Amer. Chem. Soc., 53, 3779, 1931. 
27 Emeléus, Trans. Faraday Soc., 28, 89, 1932. 
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Less than 5 per cent of the ammonia reacting yields the normal prod- 
ucts V,and H,. Since the main reaction product is deposited in con- 
densed form, the atmospheres of Jupiter and Saturn should be de- 
pleted of carbon monoxide by this time. This conclusion is the more 
remarkable since the presence of carbon monoxide is difficult to ex- 
clude on spectroscopic evidence and its vapor tension, nearly as high 
as that of nitrogen, would favor its persistence in the gaseous state. 
In addition, it may be mentioned here that the presence of carbon 
dioxide in the atmosphere of Jupiter is hard to reconcile with that of 
ammonia because these gases associate, forming ammonium car- 
baminate, 


2NH, + CO, = NH;. 


The dissociation tension of this solid compound is exceedingly small 
compared with the vapor tension of either pure ammonia or pure 
carbon dioxide. Although the dissociation tension of ammonium 
carbaminate has not yet been measured at temperatures lower than 
room temperature, it can be estimated that the surface temperature 
of Jupiter, usually assumed to be about 150° K, would have to be 
increased by at least 100° in order to maintain the spectroscopically 
observed amount of ammonia in equilibrium with the corresponding 
amount of gaseous carbon dioxide and solid ammonium carbaminate. 
The analogous argument and about the same numerical estimate ap- 
plies to the presence of hydrogen sulphide, which also would associ- 
ate with ammonia, forming VH,:SH. 


PHOTOCHEMISTRY OF SOME GASEOUS OXIDES 


The atmosphere of Venus seems to be composed almost entirely 
of carbon dioxide, and the amounts of oxygen and water vapor must 
be less than one-thousandth and one-tenth, respectively, of the 
amounts present in the earth’s atmosphere. The estimates of the 
total mass of carbon dioxide are based on the strength of the ab- 
sorption bands situated in the near infrared. They vary between 
400 meter-atmospheres and 2 mile-atmospheres. Since oxygen is 
practically absent, the carbon dioxide is exposed to the destructive 


| 
| | 
| 


PLANETARY ATMOSPHERES 335 


action of ultraviolet sunlight, which brings about the following pri- 
mary process: 


< 1800 A: CO, + hv -O+CO. 


In the highest atmospheric layers carbon monoxide may dissociate 
photochemically in its turn: 


<1200A: 


but throughout the bulk of the atmosphere carbon monoxide will be 
protected by the absorption of carbon dioxide. The most important 
secondary reaction should be the partial recombination of the oxygen 
atoms, which would be followed by a speedy photochemical dissocia- 
tion of the oxygen molecules. Carbon monoxide molecules are not 
known to associate in any way. Thus the problem is reduced to the 
oxidation of carbon monoxide by photochemically produced oxygen 
atoms, which has been studied by several investigators. The exten- 
sive work by Popov’ has little bearing on the astrophysical problem 
because in his experiments the oxidation took place under very low 
pressure, and, therefore, the homogeneous reaction was negligible 
beside that on the wall. Jackson”? used pressures up to several hun- 
dred millimeters, and his experiments cover a wide range of tempera- 
tures (300°-700° K). At room temperature oxygen atoms react 
with oxygen molecules, forming ozone at least one hundred fifty 
times more rapidly than with carbon monoxide molecules, forming 
carbon dioxide. At high temperatures the dry oxidation is practical- 
ly independent of temperature and total pressure. With water vapor 
present, the reaction is no longer independent of pressure and has a 
large temperature coefficient. A mechanism involving OH and H 
was suggested by Jackson in order to explain the effect produced by 
water vapor. Recently it has been proved* that the photochemical 
oxidation of carbon monoxide is a chain reaction, and various reac- 
tion schemes have been suggested. A final choice between these is 
not yet possible. Still it may be safely concluded that the photolysis 


28 Acta physicochemica U.S.S.R., 3, 223, 1935- 
29 W. F. Jackson, J. Amer. Chem. Soc., 56, 2631, 1934. 
3° L. Holzapfel, Inaugural-Dissertation, Universitat Berlin, 1936. 
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of carbon dioxide in the atmosphere of Venus is followed by a re- 
oxidation of the carbon monoxide. The presumably high tempera- 
ture at the surface of Venus—on account of the greenhouse effect of 
the carbon dioxide—and the presence of traces of water vapor, which 
cannot be ruled out at present, may be helpful agents. Water vapor 
itself decomposes photochemically according to the following equa- 
tion: 


< 1800A: 


Since the limiting wave-length of photolysis is approximately the 
same for water vapor and carbon dioxide, these compounds should 
be decomposed simultaneously in an atmosphere containing both of 
them. This may lead to a subsequent formation of formaldehyde, 
which in laboratory experiments has been observed to originate from 
the reaction between carbon monoxide and photochemically pro- 
duced hydrogen atoms. The same process may occur in the atmos- 
phere of Venus. Formaldehyde has been detected by chemical meth- 
ods in the earth’s atmosphere,*" but its photochemical origin has not 
yet been ascertained there. Formaldehyde displays a strong elec- 
tronic absorption spectrum in the astronomically accessible ultra- 
violet. Therefore it would be of interest to look for it in the spectrum 
of Venus. 

Laboratory experiments give some evidence of the photochemical 
oxidation of nitrogen. Oxidation of nitrogen is a large-scale phe- 
nomenon in terrestrial thunderstorms. But apparently the oxides 
are destroyed rapidly, since their stationary concentration is almost 
immeasurably small. Under these circumstances it is improbable 
that nitrogen oxides should be detected spectroscopically in the at- 
mospheres of the other terrestrial planets, all of which contain traces 
of oxygen, at best. 


PRINCETON UNIVERSITY OBSERVATORY 
June 1937 


3tN. R. Dhar and A. Ram, Nature, 132, 819, 1933. 
32,V. Kondratjev, Acta physicochemica U.S.S.R., 3, 247, 1935. 
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NOTE ON THE TEMPERATURE OF VENUS 
ARTHUR ADEL 


ABSTRACT 

The distribution of intensity of absorption in the fine structure of the CO. bands in 
the spectrum of Venus provides an independent approach to the problem of the planet’s 
temperature. The temperature of the surface appears to be greater than 50° C. 

Information regarding the temperatures of planetary atmospheres 
is generally sought through the determination of the distribution of 
energy in the reflection and emission spectra of the planets. There is, 
however, another method which is not so incumbered by the effects 
of the earth’s atmosphere and which is quite ideally suited to those 
planets whose spectra include absorption bands characteristic of 
their atmospheres: Venus contains carbon dioxide; Jupiter and 
Saturn, methane and ammonia; Uranus and Neptune, methane. 
For, in a planetary atmosphere, the distribution of molecular popu- 
lation over the various vibrational and rotational states is strictly 
determined by the prevailing temperatures. In consequence the dis- 
tribution of intensity of absorption in the fine structure of the bands 
is also governed by the temperature, and is in most instances a very 
simple function of the latter. 

*For the solution of the problem one requires, however, definitive, 
high-dispersion spectra to make possible the unambiguous selection 
of the fine-structure lines displaying the maximum intensity of ab- 
sorption. With the possible exception of a reproduction of a portion 
of the spectrum of Venus, such data are not available.’ The spec- 
trum alluded to, while scarcely adequate for the present need, will 
suffice for purposes of illustration. 

The reproduced region of the spectrum of Venus displays the two 
carbon dioxide absorption bands 


*W. S. Adams and T. Dunham, Jr., Pub. A.S.P., 44, 243, 1932. 
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and 


with centers at 12672.4 and 12774.7 cm™, respectively.? Each of 
these bands possesses a converging high-frequency branch and a 
diverging low-frequency branch, but no zero branch. Apropos of the 
congestion in the high-frequency branches, the diverging branches 
are more convenient for study. Moreover, the short wave-length 
band (12774.7 cm™) is apparently the more easily examined of the 
two. The theoretical intensity of the jth line in the negative branch 
of this band is given by the expression 


[= fou = , 


where 

ART’ 
A being the moment of inertia of the carbon dioxide molecule. The 
theoretical intensity cannot be directly observed, but the minimum 
percentage transmission can be.’ The minimum percentage trans- 
mission 7\,;, and the theoretical intensity are related by the expres- 
sion 
d{log Tmin} = , 


which applies to a layer of the Venusian atmosphere of thickness dx 


and temperature T,. The value of 7, determines the value of 7 for 
which the function Ti, is a minimum. For the entire atmosphere 


we have 


{log Timin}? = e , 


It is sufficient to point out that this function (or T,,;, for the entire 
atmosphere) will also be a minimum for a certain value of 7 and that 


2 Arthur Adel and D. M. Dennison, Phys. Rev., 43, 7163 44, 99, 1933. 
3D. M. Dennison, ibid., 31, 503, 1928. 
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this will correspond to a value of T lower than the surface tempera- 
ture of the planet and higher than the temperature of the outer lay- 
ers of the planet’s atmosphere. In other words, we may write 


{log Vital? = Ci 


where o, contains this effective temperature. 

From the water-cell transmissions made in 1922 by W. W. Cob- 
lentz and C. O. Lampland,* D. H. Menzel’ has calculated an effec- 
tive noonday temperature of 50° C. for the surface of Venus. We 
should therefore expect to find the value of 7 corresponding to 50° C. 
greater than the value of 7 at which maximum absorption occurs in 
the negative branch of 


From the expression 


it follows that the maximum intensity of absorption relates the effec- 
tive temperature with the ordinal number of the line in question in 
the manner 


= VART . 
For T = 273° + 50° = 323°,,andA = 70.2 X 10 “gm cm’, we find 
j & 18. Since alternate fine-structure lines are completely forbidden 
by the selection rules, we expect to find the maximum intensity of 
absorption somewhere between the center of the band and the ninth 
appearing line. Although it is difficult to ascertain from the repro- 
ductions exactly where the maximum absorption occurs, careful ex- 
amination indicates that 50° C. is probably too low a value for the 
temperature of the surface of the planet Venus. 


LOWELL OBSERVATORY 
FLAGSTAFF, ARIZONA 


July 1937 


4 Pop. Astr., 30, 551, 1922. SAp. J., 58, 65, 1923. 
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THE NUCLEAR STAR CLUSTER IN 30 DORADUS 


Recent photographs made with the 60-inch reflector of the Boyden 
Station permit a further description of the gaseous nebula 30 Dora- 
dus, which shares with the ‘“‘axis” the distinction of being the most 
conspicuous feature of the Large Magellanic Cloud. In particular 
the photographs reveal the nature and richness of the remarkable 
cluster of blue stars that lies in the heart of the nebulosity. In neither 
linear dimensions nor total luminosity is 30 Doradus closely ap- 
proached by any gaseous nebula of the galactic system. The nebu- 
losity, NGC 5195, at the end of a spiral arm of Messier 51, may be, 
however, of comparable dimensions, luminosity, and general struc- 
ture. 

The position for 1900.0 of 30 Doradus, which is known also as the 
looped nebula and as NGC 2070, is a= 5394, 6= —69°9’, A= 246°, 

=—31°. 

In spectral characteristics 30 Doradus resembles the great Orion 
nebula (NGC 1976). The nebular lines, N, and N,, and the hydrogen 
lines predominate, but a close study of the spectrum as a whole is 
not possible because of the large dimensions of the nebula. An exam- 
ination by Miss Cannon shows faintly and diffusely the O 1 pair of 
lines at \ 3727, not hitherto recorded for 30 Doradus, but conspicu- 
ous in the Orion nebula. 

The nebulosity, both dark and bright, extends to great distances 
from the center. The total diameter is about 25’, corresponding to 
600 light-years, which is about one hundred times the diameter of 
the Orion nebula and perhaps five times the diameter of W. H. 
Pickering’s spiral of faint nebulosity that spreads throughout the 
whole constellation of Orion. The variable HV 2740, which is 5’ from 
the center, is apparently irregular; an earlier indication that it was a 
Cepheid with period of 23.0 days, showing a large deviation from the 
period-luminosity curve," has been rejected after further study. Half 


* Harvard Bull., No. 883, 1931. 
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PLATE VII 


South 


30 DoraApbus, IN THE LARGE MAGELLANIC CLOUD 


Photographed by J. S. Paraskevopoulos, with the 60-inch reflector at the Boyden Station of Harvard Observatory. 
nlarged 2.7 times. Scale, 1 mm=70"2. Exposure, 50 minutes. 
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a dozen other variables have been found within the bounds of the 
nebula by Mrs. S. F. Lindsay, but types have not yet been deter- 
mined. The total photographic absorption produced by the nebula 
probably exceeds one magnitude. 

Five clusters of the galactic type are apparently involved in the 
nebulosity. The richest is the loose, elongated cluster at the center. 
Its major and minor axes are approximately sixty and forty-five 
light-years. Estimates of magnitude are difficult because of the ir- 
regular nebulous background and the Eberhard effect. An examina- 
tion of the nuclear cluster on seven plates of different exposure-times 
gives for the number of stars brighter than various magnitude limits: 


A comparison of photographs on blue- and red-sensitive plates 
indicates that all of the cluster stars are blue. Their exceedingly high 
luminosities are, of course, responsible for the incited radiation of 
the large nebula in which they are involved. Even with no allow- 
ance made for the strong absorption within the nebula, their abso- 
lute magnitudes range from positive values up to —5 and brighter. 

In the nuclear area Miss Cannon has classified the spectra of five 
stars.? Some are fainter than the fourteenth magnitude, but all are 
of spectral class O. 

It is of interest to note that, notwithstanding the large number 
of stars in the nuclear cluster and their exceedingly high luminosity, 
they contribute altogether but a small portion—perhaps 10 per cent 
—of the total light of 30 Doradus, which is a naked-eye object. 
Estimates made at the Lick Observatory’ and at Harvard place the 
integrated apparent magnitude of the nucleus stars and nebulosity 
at about 6.0; a Harvard estimate of the integrated magnitude of the 
whole nebula, with the small contribution from the involved clusters, 
gives the value 4.0.4 > 

2 Ibid., No. 891, 1933; Harvard Ann., 76, 28, 1916. 

3 Lick Obs. Pub., 13, 109, 1918; see also ibid., pp. 187-80. 


4 Harvard Circ., No. 271, 1925. In this paper, through a misprint, the total diameter 
is given as 4o parsecs; it should read 140 parsecs, the value indicated by the plates then 
available. 
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The accompanying photograph was made on January 17, 1937, 
with an exposure time of fifty minutes on a Cramer Hi-Speed Special 
plate. The original photograph shows objects fainter than the 


nineteenth magnitude. 
HARLOW SHAPLEY 
JOHN S. PARASKEVOPOULOS 


HARVARD OBSERVATORY 
September 1937 


THE PHOTOELECTRIC COLOR OF THE 
ZODIACAL LIGHT 


ABSTRACT 


The mean color of the zodiacal light as observed with a photoelectric photometer at 
McDonald Observatory is +0.28 mag. on Becker’s system of photoelectric colors. This 
color corresponds to spectral type G1, which is very near the color of the sun. It is con- 
cluded that there is no Rayleigh scattering in the zodiacal cloud and hence the number 


of very small meteoric particles is negligible. 

Colorimetric measures of the zodiacal light are data which will 
be useful in an analytical discussion of the physical properties of the 
medium producing the illumination. It is generally believed that the 
zodiacal light is produced by the scattering of sunlight by meteoric 
particles concentrated in the plane of the ecliptic and extending be- 
yond the orbit of the earth. According to the theories of scattering, 
if the particles are small—say one-quarter of the wave-length of 
light—the radiation scattered by them will follow Rayleigh’s law, 
that is, the zodiacal light in this case should be as much bluer than 
the sun as is the daytime sky. If the particles are large compared 
with the wave-length of light, then the scattering is that of pure re- 
flection and is independent of the wave-length of the incident light. 
For intermediate sizes the law of scattering is somewhat complicated 
and is a function of various powers of the wave-length. 

The color of the zodiacal light has been measured by one of us’ 
with a photoelectric photometer at the Yerkes Observatory. The 
results were not consistent, indicating that possibly the color of the 
zodiacal light changed from time to time. It was stated, however, 
that a part of this variation might be due to the light from the night 
sky, since the colorimetric determinations were for the total light, 
that is, the auroral light, the scattered light, and the light of the un- 
resolved stars, plus the zodiacal light. 

Ap. J., 80, 61, 1934. 
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In view of the large variations found in that work it seemed wise 
to take additional observations of the color of the zodiacal light un- 
der better conditions for measurements at such large zenith dis- 
tances. We have observed this color at the McDonald Observatory, 
with the photoelectric photometer used and described by Elvey and 
Roach.? Owing to the high elevation, 2070 m, the sky is quite uni- 
form and free from haze at large zenith distances. The method of 
observing was quite different since with the alt-azimuth arrangement 
of the photoelectric sky photometer it was quite easy to make a se- 
ries of observations across the zodiacal-light cone at a given zenith 
distance. The observations near the ends of the run gave the sky 
correction to be applied to the zodiacal light, since the auroral and 
scattered light is constant for a given zenith distance. In this man- 
ner nothing but the zodiacal light was included in the observations 
for color. Since the observations were taken sufficiently far from the 
Milky Way, there was very little differential effect from the unre- 
solved stars between the areas used for the sky background and the 
region of the zodiacal light. The colors were calibrated and cor- 
rected for extinction by observations on eight bright stars of various 
spectral classes. This number of stars is not sufficient to establish 
independently the relationship between color and spectral type, but 
there are enough data to determine the color equation between our 
observed colors and the photoelectric colors of stars measured by 
W. Becker. The resulting color equation for reducing our observa- 
tions to Becker’s system is 


Color (Becker) = 0.37 + 0.74 X Color (obs) . 


The accompanying tabulated results are for the evening zodiacal 
light. The mean color of the zodiacal light is —o.12 mag, which 


Date Color 


2 Tbid., 85, 213, 1937. 
3 Veriff. Sternwarte Berlin-Babelsberg, 10, Nos. 3 and 6, 1933. 
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when reduced to Becker’s system is equal to +0.28 mag. This color, 
! when referred to his diagram showing the dependence of color upon 
spectral type, corresponds to spectral class Gr. 

Since the foregoing color of the zodiacal light is very similar to 
that of the sun, we can conclude that there is practically no Ray- 
leigh scattering in the cloud of meteoric material and hence the num- 


ber of very small particles is negligible. 
C. T. ELVEY 
PAUL RUDNICK 
McDoNALpD OBSERVATORY 
Fort Davis, TEXAS 
July 23, 1937 
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